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ABSTRACT 


An  analysis  of  numerical  methods  for  extracting 
aerodynamic  coefficients  from  dynamic  test  data  has  been 
conducted.  The  emphasis  of  the  analysis  is  on  the  effects 
that  random  measurement  errors  in  the  data  and  random 
disturbances  in  the  system  have  on  the  accuracy  with  which 
the  coefficients  for  linear  and  nonlinear  systems  can  be 
determined.  Both  deterministic  and  stochastic  methods  for 
extracting  the  coefficients  and  determining  their  uncer¬ 
tainties  are  considered. 

The  deterministic  technique  considered,  due  to 
Chapman  and  Kirk,  provides  excellent  estimates  of  both 
linear  and  nonlinear  static  pitching  moment  coefficients 
for  the  range  of  measurement  errors  and  system  noise 
considered.  Somewhat  less  accurate  estimates  of  linear 
damping  coefficients  are  obtained.  Nonlinear  pitch  damp¬ 
ing  coefficients  extracted  using  this  deterministic  technique 
are  affected  considerably  by  both  measurement  errors  and 
system  noise.  The  estimated  standard  deviations  of  the 
extracted  coefficients  obtained  using  standard  techniques 
are  generally  adequate  when  the  data  being  analyzed  con¬ 
tain  only  measurement  errors. 

The  stochastic  approach  considered  demonstrates  the 
feasibility  of  using  an  extended  Kalman  filter,  with  a 
parameter  augmented  state  vector,  for  determining  the 
values  of  the  aerodynamic  coefficients  and  their  uncer¬ 
tainties  from  noisy  dynamic  test  data. 

The  specific  filter  used  generally  reaches  near 
steady-state  conditions  in  its  estimates  of  the  parameters 
in  less  than  one  second.  Variations  in  the  initial  para¬ 
meter  variances  or  in  the  estimates  of  the  noise  statistics 
essentially  affect  only  the  determination  of  the  nonlinear 
damping  parameter. 

Parameter  estimates  obtained  from  the  extended  filter 
compare  favorably  with  previously  obtained  results  using 
deterministic  techniques .  Estimates  of  the  parameter 
uncertainties  provided  by  the  filter  are  generally  superior 
to  those  obtained  with  deterministic  techniques  particularly 
when  system  noise  has  corrupted  the  data. 
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SECTION  I 


INTRODUCTION 


1.  Subject  Matter 

This  report  presents  an  analysis  of  methods  for 
determining  aerodynamic  coefficients  of  flight  vehicles 
for  dynamic  test  data.  This  determination  is  accomplished 
by  finding  numerical  values  of  the  aerodynamic  coefficients 
appearing  in  the  equations  of  motion  such  that  the  solu¬ 
tions  to  these  equations  are  adequate  representations  of 
the  given  test  data.  When  this  determination  has  been 
made,  the  coefficients  are  said  to  have  been  extracted 
from  the  dynamic  data. 

The  emphasis  of  the  analysis  is  placed  on  the  effect 
that  random  measurement  error-  in  the  data  and  random 
disturbances  in  the  system  have  on  the  accuracy  with  which 
the  coefficients  for  linear  and  nonline  -r  systems  can  be 
determined.  Both  deterministic  and  stochastic  system 
models  for  extracting  the  coefficients  and  determining 
their  uncertainties  are  considered. 

2.  Historical  Background 

The  determination  of  aerodynamic  coefficients  from 
dynamic  data  is  generally  agreed  to  have  begun  with  the 
work  of  Fowler,  Gallop,  Lock,  and  Richmond 1  in  the  first 
quarter  of  this  century.  Their  basic  technique  was  con¬ 
cerned  with  determining  moment  and  damping  characteristics 
of  artillery  shells  by  firing  these  projectiles  through 
spaced  cards  and  reconstructing  the  pitch  and  yaw  angle 
time  histories  by  observing  the  obliqueness  of  the  holes 
that  resulted  when  the  shells  passed  through  the  cards. 
Their  technique  of  data  measurement  is  still  in  use  at 
some  ballistic  ranges  to  this  date.*2 

.  3 

Nielsen  and  Synge  later  clarified  the  linear  theory 
of  Fowler  et  al  in  their  work  during,  and  immediately 
following,  World  War  II. 

Coefficient  extraction  techniques  that  are  currently 
in  use  at  many  ballistic  and  wind  tunnel  facilities  evolved 
from  the  work  of  Murphy 4, 5  and  Nicolaides,6  both  of  whom 
have  considered  various  combinations  of  degrees  of  freedom 
for  a  variety  of  flight  vehicles.  Their  aerodynamic  co¬ 
efficient  extraction  techniques  feature  a  least  squares 
fit  to  dynamic  data  of  the  exact  solutions  of  linear 
equations  of  motion,  or  approximate  closed  form  solutions 
of  slightly  nonlinear  equations  of  motion  using  the  quasi¬ 
linearization  technique  of  Kryloff  and  Bogoliuboff J  Ex- 
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tensions  of  the  work  of  Murphy  and  Nicolaides,  particularly 
for  more  complex  nonlinear  systems,  have  been  made  by 
Eikenberry  and  Ingram. 9 

The  requirement  for  closed  form  solutions  of  the 
equations  of  motion  has  recently  been  eliminated  by  the 
formulation  of  least  squares  techniques  which  fit  numerical 
solutions  of  the  equations  of  motion  to  dynamic  data.  The 
minimization  of  the  least  squares  criterion  function  involves 
differential  corrections,  as  do  the  techniques  of  Murphy 
and  Nicolaides;  the  required  partial  derivatives  are 
determined  by  numerically  integrating  parametric  differ¬ 
ential  equations  which  are  derived  from  the  equations  of 
motion.  Contributions  to  this  numerical  coefficient 
extraction  technique  have  been  made  independently  by 
Chapman  and  Kirk,  10  Knadler,H  Goodman^,  13  and 
Meissinger ,14  although  the  method  is  usually  referred  to 
as  the  Chapman-Kirk  technique.  The  computational  require¬ 
ments  of  this  technique  are  sometimes  extensive, 15  but 
this  is  usually  outweighed  by  the  fact  that  it  can  be 
used  to  analyze  highly  nonlinear  aerodynamic  forces  and 
moments . 16 

All  of  the  coefficient  extraction  techniques  that 
have  been  discussed  to  this  point  are  deterministic  in 
nature  in  that  the  modeling  of  the  equations  of  motion 
does  not  account  for  random  disturbances  in  the  system. 

Most  angular  or  translational  motion  data  obtained  from 
ballistic  ranges  or  wind  tunnel  dynamic  tests  have,  never¬ 
theless,  been  affected  by  random  system  noise  and  random 
measurement  errors.  The  effects  of  noise  of  these  types 
on  the  accuracy  of  coefficients  extracted  from  dynamic  data 
using  deterministic  techniques  are  of  current  interest. 17 
With  the  exception  of  Ingram's^  partial  analysis  of  noise 
effects  on  Eikenberry' s  "Wobble"  program, 8  however,  little 
has  been  reported  concerning  these  effects.  One  of  the 
primary  goals  of  the  research  reported  here  is  to  show  some 
effects  of  random  measurement  errors  and  system  noise  on 
the  accuracy  of  aerodynamic  coefficients  extracted  from 
dynamic  data  using  the  most  general  of  the  various  determi¬ 
nistic  techniques,  ’:.hat  of  Chapman  and  Kirk. 

In  recent  years  increasing  attention  has  been  given 
to  parameter  and  state  variable  estimation  through  the 
use  of  stochastic  modeling  of  the  physical  system  of 
interest.  Most  of  this  work  has  been  done  by  optimal 
control  specialists  and  is  a  direct  result  of  pioneering 
contributions  in  linear  filtering  theory  by  Kalmanl8,19 
and  Bucy.19  The  Kalman  filter  provides  estimates  of  the 
states  of  noisy,  linear  dynamic  systems  as  well  as  estimates 
of  the  state  variable  uncertainties.  Bryson  and  Ho, 20 
among  others,  have  extended  the  Kalman  filter  theory  to 
include  estimates  of  the  states  for  nonlinear  systems. 
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Mehro^l  has  recently  proposed  a  maximum  likelihood 
technique  for  the  determination  of  aerodynamic  coefficients 
from  dynamic  data  using  the  Kalman  filter  fox  linear  systems 
and  the  extended  Kalman  filter  when  the  system  is  non¬ 
linear.  Results  verifying  his  proposals  are,  as  yet, 
unavailable. 

An  additional  maximum  likelihood  technique  has  been 
developed  by  Grove  et  al22  and  has  been  used  with  modest 
success  by  Suit. 23 

A  stochastic  approach  to  the  coefficient  exc. 'action 
problem  using  an  extended  Kalman  filter  with  paranu.  ser 
augmented  state  vector  is  developed  in  this  report. 

Results  obtained  with  this  filter  are  also  presented. 

3 .  Scope  of  the  Report 

The  remainder  of  this  report  is  devoted  to  an  analysis 
of  the  effects  of  random  system  noise  and  measurement 
errors  on  coefficients  extracted  from  dynamic  data  using 
the  method  of  Chapman  and  Kirk,  as  well  as  the  development 
and  evaluation  of  an  extended  Kalman  filter  for  solving 
essentially  the  same  problem.  The  equation  of  motion 
used  in  the  analysis  is  one  describing  che  one-degree- 
of-freedom  pitching  motion  of  a  rigid  body.  Pitching 
moments  that  are  both  linear  and  nonlinear  functions  of 
angle-of-attack  are  considered. 

The  theory  of  the  Chapman-Kirk  technique  is  recounted 
briefly  in  Section  II,  followed  by  the  analysis  of  noise 
effects  on  this  method  in  Section  III. 

Developments  leading  to  the  formulation  of  the  extended 
Kalman  filter  to  be  used  are  given  in  Section  IV.  Results 
obtained  with  the  extended  filter  are  presented  and  discussed 
in  Section  V. 

A  summary  is  given  in  Section  VI. 

Computer  program  listings  and  definitions  of  program 
variables  are  given  in  the  two  appendices. 


SECTION  II 


THE  COEFFICIENT  EXTRACTION  TECHNIQUE 
OF  CHAPMAN  AND  KIRK 


Development  of  the  Extraction  Algorithm 


A  brief  reconstruction  of  the  Chapman-Kirk  algorithm 
for  extracting  aerodynamic  coefficients  from  dynamic  data 
is  given  in  this  section.  The  standard  technique  for 
estimating  uncertainties  in  coefficients  determined  from 
a  least  squares  fit  of  a  given  function  to  experimental 
data  is  also  presented. 

The  extraction  technique  of  Chapman  and  Kirk  has  two 
very  basic  requirements:  (1)  the  differential  equation 
of  motion  for  the  body  of  interest  must  be  given  and 
(2)  a  set  of  experimental  data  based  on  the  observed  motion 
of  the  bcdy  must  be  available. 

As  an  example,  consider  the  nonlinear  equation  of 
pitching  motion  for  a  rigid  body 

fcl  +  (C4  +  C0ct2)a  +  (C3  +  Csa2)a  =  0 

subject  to  the  initial  conditions 

a(0)  =  Ci 

a (0)  =  C2 

where  a  indicates  the  derivative  of  a  with  respect  to  time. 
Suppose,  also,  that  the  time  history  of  the  pitch  angle 
as  recorded  during  some  experiment,  a  (t^) ,  i=l,2,...m, 
is  also  available.  e 

The  technique  of  Chapman  and  Kirk  is  used  to  determine 
the  values  of  the  Cj  (j=l,2, . . . ,6)  in  Equation  (1)  which 
result  in  the  solution  to  this  equation  of  motion  be? ng 
a  best  fit  to  the  test  data  in  a  least  squares  sense. 

Thus,  it  is  necessary  to  minimize  the  least  squares  cri¬ 
terion  function 


=ii  -°e<ti>  '  VV 


where  a  (t. )  is  obtained  from  the  solution  to  Equation  (1) . 

C  i. 
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Nov?  it  is  well  known  that  in  order  to  determine  para¬ 
meters  directly  by  a  least  squares  fit  of  a  given  function 
to  test  data,  the  parameters  must  appear  linearly  in  the 
function.  This  requirement  is  met  in  the  problem  at  hand 
by  expanding  ac  in  a  truncated  Taylor  series  about  the 
numerical  solution  resulting  from  some  initial  estimates 
of  the  parameters  of  interest,  C.  i.e., 


ac(ti> 


a  (t . 
co  1 


i!  +^fe)i 


where 


(  3ct 

\  ^ 


\  are  evaluated  for  C.=C.  and  AC.=C.-C.  . 
j  3  Do  3  3  30 


Substituting  (3)  into  (2)  yields 


s  -  J  ..<**> 


aco(ti^  3C . 

3  3 


Now,  assuming  that  the 


3a 

3C. 

L  3  J 


are  known.  Equation  (4)  is  a 


function  of  the  AC j 1 s  only ,  Therefore ,  to  determine  the 
values  of  these  coefficients  that  will  minimize  this  equation, 
it  is  necessary  to  take  the  partial  derivative  of  Equation 
(4)  with  respect  to  each.  AC.,  set  each  of  the  resulting 
equations  to  zero  and  solve-5 for  the  ACj's.  Carrying  out 
these  operations  yields  the  following  matrix  equation 


[A]  [AC]  =  [B] 


where  A  is,  in  this  case,  a  6  X  6  matrix  with  elements 
given  by 


n  r 

a  _  3a  3a 

jk  =  iii  L  acj  J  4  L3Ck  J  i  ' 
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AC  is  a  6  X  1  column  matrix,  or  vector,  and  B  is  a  6  X  1 
column  matrix  with  elements  given  by 


B . 
3 


m 

I 

i-1 


ae(ti} 


“co^i* 


( 


The  solution  to  Equation  (5)  is 


[AC] 


[A]-1 


[B] 


( 


The  solutions  for  the  ACjls  obtained  from  Equation  (8) 
are  exactly  correct  only  if  ac  is  a  linear  function  of  the 
Cj's  as  assumed  by  Equation  (3).  This  condition  of  linearity 
is  seldom  the  case,  and  the  process  must  be  repeated  with 
new  initial  guesses , 


C.  +  AC. 
DO  D  o 


until  the  change  in  the  criterion  function  [Equation  (2) ] 
from  one  iteration  to  the  next  is  sufficiently  small. 

The  algorithm  just  presented  requires  that  time 
histories  of  the  influence  coefficients,  3g  ,  be  available. 

dC . 

D 

These  time  histories  are  determined  by  numerically  inte¬ 
grating  parametric  differential  equations  which  are  derived 
by  differentiating  the  equation  of  motion  with  respect  to 
each  of  the  parameters  of  interest.  As  an  example,  the 
parametric  differential  equation  for  Ci,  the  initial 
condition  of  a  is 


3CX 


[8  +  (C4+C6a2)a  +  (C3+C5a2)ct]  =  0 


Itv-  +  ivr*-  a2a  +  2C6aa  ~r~  +  C6cx2  t^—  +  a 


3Ci  3CX 


9C  i 


3Cx 


acx  dCi 


+  C3  +  Ir5-  a2  +  3C5a2  =  0 


T 


Assuming  that  the  parameters  are  independent  of  each  other 
and  that  the  order  of  differentiation  can  be.  reversed, 
the  final  form  of  the  desired  equation  is 


dt2  I  a Cj 


d2  3a  ,  -  o  d 
-  -  +  c6a-  - 


— —  — —  +  (C3+3C5a2+2Cgaa)  ^2- 

dt  LsclJ  LSCi 


subject  to  the  initial  conditions 


3a  (0)  _  3Cj 


3a  (0) 


The  complete  set  of  parametric  differential  equations  for 
the  given  equation  of  motion  [Equation  (1)]  is  given  in 
Section  III. 

In  summary,  the  process  for  extracting  numerical 
coefficients  from  test  data  given  the  system  model  [Equa¬ 
tion  (1)]  and  criterion  function  [Equation  (2)]  is 

1.  Estimate  the  numerical  values  of  the  Cj's. 

2.  Integrate  Equation  (1)  to  obtain  acofti) • 


3 .  Determine 


rial  . 
UcjJi 


4.  Solve  Equation  (8)  for  the  ACj's. 

5.  Repeat  the  process  with  Cj i=Cjp=ACj0  until  the 
change  in  Eugation  (2)  is  sufficiently  small. 
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2. 


Estimation  of  Extracted  Parameter  Uncertainties 


The  estimation  of  the  uncertainties,  or  standard 
deviations,  of  parameters  that  have  been  determined  by  the 
lec  _  squares  fitting  of  a  given  function  to  test  data  is 
a  well-known  result  available  in  a  variety  of  references 
(see,  for  example,  References  24  or  25) . 

The  least  square  parameter  uncertainties  are  estimated 
in  this  report  by 


a 


j 


(12) 


where  Aj-j-1  is  the  jth  diagonal  element  of  the  inverse 
Grammian  matrix  [Equations  (6)  and  (8)],  S  is  the  sun.  of 
the  squares  of  the  residuals  as  given  by  Equation  (2) ,  m 
is  the  total  number  of  data  points,  and  K  is  the  total  number 
of  parameters  being  determined  by  the  fit. 
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SECTION  III 


ANALYSIS  OF  THE  CHAPMAN-KIRK 
COEFFICIENT  EXTRACTION  TECHNIQUE 

1.  Coefficient  Extraction  Computer  Program 

This  section  provides  a  detailed  description  of  a 
one-degree-of-freedom  coefficient  extraction  computer 
program  based  on  the  previously  described  iterative  process 
of  Chapman  and  Kirk.  This  program  considers  the  pitching 
motion  of  a  symmetric  missile  about  a  fixed  point.  It 
is  used  to  determine  values  of  static  pitching  moment 
coefficient  derivatives,  pitch  damping  coefficients,  and 
a  trim  term  so  that  the  solution  to  the  appropriate  differ¬ 
ential  equation  of  motion  is  a  best  fit  to  test  data  in  a 
least  squares  sense. 

The  program  was  developed  to  be  used  as  an  economical 
tool  in  the  evaluation  of  the  sensitivity  to  noise  and 
convergency  sensitivity  of  the  Chapman-Kirk  technique  when 
operating  in  its  least  complex  mode. 

a.  Computational  Equations 

(1)  Equation  of  Motion 

The  complete  equation  of  motion  that  is  used  in  the 
program  is 


a  +  (C4+C6a2)a  +  (C3+C5a2+C7a4)a  +  C8  =  0 


(13) 


with  initial  conditions 


o(0)  -  Ci  ,  d(0)  =  C2 


(14) 


where 


C3 


(C  ) qAd 
mao  ’ 

I 


C4 


(C  ) qAd2 
mqo  ^ 

2VI 


9 


(Cma2)9M 


(cmc.4)<aA<a 

I 


(Cmq2)qM2 


(Cm{g)qad 

I 


(2)  Parametric  Differential  Equations 

The  eight  parametric  differential  equations  for  the 
given  equation  of  motion  [Equation  (13) ]  and  initial  con¬ 
ditions  [Equation  (14)]  are  of  the  form 


P.  +  AP.  +  BP.  =  F.  ,  .  =  1,2,..., 8 


where 


p  -  d_ 
j  dt  3Cj 


-  (C4  +  Cgct2) 


p  _  d2  3a 
j  dt7"  3Cj 


B  =  (C3  +  3C5a2  +  2C6aa  +  5C7a4) . 


The  values  of  the  initial  conditions  Pj  (0)  and  Pj  (0)  as 
well  as  the  functional  forms  of  the  nonhomogeneous  term 
Fj  are  given  in  Table  I  for  j=l,2,...,8. 
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TABLE  I.  INITIAL  CONDITIONS  AND  NONHOMOGENEOUS 

TERMS  FOR  PARAMETRIC  DIFFERENTIAL  EQUATIONS 


j 

Pj  (0) 

Pj(0) 

F. 

3 

1 

1 

0 

0 

2 

0 

1 

0 

1  3 

0 

0 

-a 

4 

0 

0 

• 

-a 

5 

0 

0 

-a3 

6 

0 

0 

-aa2 

7 

0 

0 

i 

a 

8 

0 

0 

-1 

b.  Program  Description 

The  program  is  written  in  Fortran  IV  for  use  primarily 
on  an  IBM  360/65.  The  paragraphs  that  follow  describe  the 
functions  of  the  main  program  and  its  four  subroutines , 
the  required  imput  data,  and  the  program  options.  A  listing 
of  the  complete  program  is  given  in  Appendix  I . 

(1)  Main  Program  Functions 

The  functions  of  the  main  program  in  their  approximate 
order  of  occurrence  are  as  follows: 

(a)  Reads  and  writes  input  data. 

(b)  Computes  initial  parameter  values,  [These  are  the 
Cj's  that  appear  in  Equation  (13).] 

(c)  Calls  the  numerical  integration  subroutine  ADDUM. 

(d)  Writes  current  parameter  values. 
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(e)  Computes  the  matrix  elements  required  for  incre¬ 
menting  the  parameter  values. 

(f)  Calls  the  matrix  inversion  subroutine  MINV. 

(g)  Computes  the  sum  of  the  squares  of  the  residuals 
between  the  calculated  and  experimental  data 
points . 

(h)  Computes  the  root-mean-square  residual  and  root- 
mean-square  error  (Reference  24)  of  the  current 
fit  to  date. 

(i)  Computes  the  estimated  standard  deviations 
(Reference  24)  of  the  current  parameter  values. 

(j)  Tests  the  difference  between  the  current  and 
previous  values  of  the  root-mean-square  error 
to  determine  if  the  iteration  process  has 
converged , 

(k)  Computes  the  incremental  changes  for  the  para¬ 
meters  (if  convergence  has  noc  occurred) . 

(l)  Computes  the  updated  parameter  values  (if 
convergence  has  not  occurred) . 

(m)  Returns  to  (c)  (if  convergence  has  not  occurred) . 

(n)  Computes  the  values  of  the  extracted  coefficients 
from  the  current  parameter  values  (after  con¬ 
vergence  is  achieved) . 

(o)  Computes  the  estimated  standard  deviations  of 
the  extracted  aerodynamic  coefficients  (after 
convergence  is  achieved) . 

(p)  Writes  extracted  aerodynamic  coefficients , 
estimated  standard  deviations,  and  the  pitch 
angle  output  that  represents  the  final  fit  to 
the  experimental  data  (after  convergence  is 
achieved) . 

(2)  Subroutine  Functions 

The  names  and  functions  of  each  of  the  four  subroutines 
that  are  used  in  the  coefficient  extraction  program  are 
given  below. 

ADDUM. —  Subroutine  ADDUM  integrates  the  equation 
of  motion  and  parametric  differential  equations  using  a 
fourth  order  Runge-Kutta  method  for  starting  and  a  fourth 
order  Adams-Bashforth  predictor-corrector  method  for 
running.  It  calls  subroutines  XDOT  and  OUT.  This  sub¬ 
routine  is  described  in  detail  in  Reference  26. 
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XDOT . —  Subroutine  XDOT  computes  current  values  of 
first  derivatives  that  are  required  by  ADDUM. 

OUT. —  Subroutine  OUT  stores  the  results  of  the 
numerically  integrated  equation  of  motion  [Equation  (13) ] 
and  parametric  differential  equations  [Equation  (15)]. 

MINV. —  Subroutine  MINV  inverts  a  Kx  K  matrix  using 
a  standard  Gauss-Jordan  technique  and  is  described  in 
detail  in  Reference  27. 

(3)  Required  Input  Data 

The  program  reads  six  categories  of  input  data. 

These  categories  and  the  specific  data  in  each  are  de~ 
lineated  in  Appendix  I.  The  format  and  units  of  the 
entries  on  a  specific  data  card  can  be  determined  from 
the  program  listing  r  >d  nomenclature  list  provided  in 
this  appendix. 

(4)  Program  Options 

This  program  has  options  for  extracting  various 
combinations  of  aerodynamic  coefficients  from  the  given 
test  data,  in  addition  to  the  option  of  extracting  no 
coefficients  and  merely  integrating  the  equation  of  motion. 
These  options  are  controlled  by  the  numerical  value  of  the 
number  of  first  order  equations  to  be  integrated,  N,  which 
is  read  by  the  program  on  the  first  data  card.  The  various 
options  are  given  in  Table  IX. 


N 

Coefficients  to  be  Extracted 

2 

No  coefficients  are  extracted.  The  program 
integrates  the  equation  of  motion  for  the 
given  input  data  and  prints  the  results. 

8 

a 

o, 

• 

a 

o, 

Cmao 

10 

“o, 

“o, 

C 

mao, 

C 

mqo 

12 

CJ 

o. 

• 

a 

o, 

C 

mao, 

C 

mqo. 

Sna2 

14 

a 

O, 

“o. 

C 

m«o, 

c 

mqo, 

Sraa2,  Snq2 

16 

O, 

• 

a 

o. 

C  ^ 
mao. 

C 

mqo. 

C  C 

ma2,  mq2, 

Cma4 

18 

1 

dO, 

• 

a 

o, 

^itlaO , 

C 

mqo, 

C  C 

ma2,  mq2, 

Sna4,  SnSa 

2 ,  Analysis  of  the  Sensitivity  to 

Noise  of~the  cE'apman-Kirk  Technique 

The  evaluation  of  the  sensitivity  to  noise  of  the 
Chapman-Kirk  coefficient  extraction  technique  has  been 
approached  from  several  directions .  Coefficients  have 
been  extracted  from  numerous  sets  of  computer  program 
generated  data  containing  only  measurement  errors  as  well 
data  containing  system  noise  and  measurement  errors . 

In  addition  to  the  use  of  artificially  generated  data, 
some  very  preliminary  work  has  been  done  with  noisy  wind 
tunnel  data  obtained  from  actual  experimentation.  The 
design  of  the  experiment  used  to  produce  the  dynamic 
data  is  described  by  Turner  in  Reference  28  along  with  some 
preliminary  results  obtained  witn  the  one-degree-of-freedom 
program  described  earlier  in  this  section  of  the  report. 

Since  the  wind  tunnel  experimentation  is  still  in  a 
developmental  stage,  howeyer,  the-  results  presented  in  the 


14 


remainder  of  this  section  are  those  obtained  from  arti¬ 
ficially  generated  data. 

a.  Generation  of  Noisy  Dynamic  Data 

The  computer  program  UFNOISE  described  in  detail  in 
Reference  29  was  used  to  generate  the  dynamic  data  from 
which  the  aerodynamic  coefficients  were  extracted.  This 
program  simulates  the  pitching  motion  of  symmetric  missile 
oscillating  in  a  wind  stream  about  a  fixed  point,  for  any 
given  initial  pitch  angle  displacement  and  initial  pitch 
angle  rare,  oy  numerically  integrating  the  equation  of 
motion.  The  program  has  two  options  for  simulating  system 
noise;  it  considers  the  magnitude  and  direction  of  the 
freestream  velocity  vector  as  normally  distributed  random 
variables ,  with  programmer  set  mean  values  and  standard 
deviations,  to  determine  the  system  noise  perturbation 
accelerations  or  it  simply  selects  random  perturbation 
accelerations  which  have  zero  mean,  are  normally  distributed, 
and  have  a  programmer  set  standard  deviation.  Both  methods 
are  essentially  equivalent.  New  values  of  the  system 
noise  perturbations  are  randomly  selected  at  each  numerical 
integration  step.  The  program  also  has  the  option  of 
simulating  random  measurement  errors  of  the  pitch  angle 
by  super imposing  normally  distributed  random  noise  on  the 
output  of  the  numerical  integration. 

The  basic  equation  of  motion  used  to  generate  the 
noisy  dynamic  data  was  a  slightly  simplified  form  of 
Equation  (13) : 


H  +  (C4  +  C6a2)a  +  (C3  +  C5a2)  a  =  w(t) 


(16) 


subject  to 


a  (0)  =  Cj  ,  a(0)=C2 


(17) 


The  true  values  of  the  physical  and  aerodynamic  con¬ 
stants  used  in  generating  the  data  are  given  in  Table  III. 
The  various  noise  level  standard  deviations  and  mean  values 
are  discussed  in  the  following  paragraphs  to  which  they 
are  pertinent. 
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TABLE  III.  TRUE  VALUES  OF  CONSTANTS 
USED  IN  GENERATING  DATA 


Constant 

Value 

d  (Ft) 

0.333 

A  (Ft2  ) 

0.0873 

I (slug  ft2) 

0.1080 

Vm (ft/sec) 

500.0 

q( lb/ft2) 

297.0 

C  (rad-1) 

mao 

-2.00 

C  (rad-3) 

ma2 

-24.5 

C 

mao 

-60.0 

C  (rad-2) 

mz 

-163.0 

a  (rad) 
o 

0.5235 

aQ (rad/sec) 

o 

• 

o 

b.  Effects  of  Random  Measurement  Errors 

Aerodynamic  coefficients  were  extracted  from  a  total 
of  nine  sets  of  data  containing  only  measurement  errors. 
Each  data  set  was  made  up  of  201  discrete  points  which  re¬ 
sult  from  integrating  the  equation  of  motion  [Equation  (16)] 
in  increments  of  0.005  second  for  a  total  of  1.000  second. 

The  desired  mean  value  of  the  measurement  errors  for 
each  of  the  nine  data  sets  was  0.0  radian.  The  desired 
standard  deviations  of  the  random  errors  were  0.00146 
radian,  0.00582  radian,  and  0.01745  radian,  with  three 
different  sets  of  data  being  generated  at  each  noise  level. 
These  standard  deviations  are  typical  of  measurement  un¬ 
certainty  for  data  of  this  type  (Reference  9) . 
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(18) 


The  above  standard  deviations  correspond  to  percent 
noise  levels  of  1.1,  4.4,  and  13.2,  respectively,  where 
percent  noise  level  is  defined  in  Reference  30  as 

Percent  Noise  =  _ _ 3o _ _ 

Approximate  average  peak  amplitude 

with  a  being  the  standard  deviation  of  interest.  The 
approximate  average  peak  amplitude  can  be  determined  in  a 
variety  of  ways.  The  value  used  here  is  the  mean  positive 
peak  amplitude  for  the  three  cycles  that  result  when  the 
nominal  equation  of  motion  is  integrated  for  1.000  second. 

The  results  of  this  portion  of  the  analysis  are  given 
in  Table  IV  and  in  Figures  1  through  10. 

Table  IV  is  primarily  a  summary  of  the  percent  error 
in  the  extracted  aerodynamic  coefficients  together  with 
their  normalized  estimated  standard  deviations  for  the 
various  actual  measurement  errors,  0^.  The  percent  error 
in  the  extracted  coefficients  is  defined  by 

Percent  Error  (C.)  =  (19) 

C. 


where  Cj  is  the  coefficient  estimate  determined  with  the 
program  in  the  fitting  of  the  data  and  Cj  is  the  true 
value  of  the  coefficient  of  interest.  Tne  normalized 
estimated  standard  deviations  are  the  estimated  standard 
deviations  calculated  from  Equation  (12)  divided  by  the 
true  value  of  the  parameter  of  interest, 


V  • 

Percent  a.  =  ^  •  100  (20) 

in  C . 


Table  IV  also  contains  the  RMS  residual  for  each  fit  to  the 

noisy  data  and  the  percent  noise  levels  based  on  the  values 

of  a- 
m 


TABLE  IV.  SUMMARY  OF  RESULTS  FOR  EXTRACTING  AERODYNAMIC  COEFFICIENTS 
FROM  DYNAMIC  DATA  CONTAINING  RANDOM  MEASUREMENT  ERRORS 


An  example  of  one  of  the  program  fits  to  a  noisy  data 
set  is  shown  on  Figure  1.  Figures  2  through  5  show  the 
variation  of  the  percent  error  in  the  individual  extracted 
coefficients  wirh  measurement  error.  The  variation  of  RMS 
residual  with  measurement  error  is  shown  on  Figure  6 .  The 
variations  of  the  normalized  estimated  standard  deviations 
of  the  extracted  coefficients  with  measurement  error  are 
shown  in  Figures  7  through  10. 

An  analysis  of  Table  IV  and  Figures  2  through  10 
reveals  the  following: 

a.  The  extracted  static  pitching  moment  coefficient 
derivatives  show  little  or  no  error  for  the 
entire  range  of  measurement  errors  considered. 

(See  Figures  2  and  3.) 

b.  The  extracted  pitch  damping  coefficients  show 
some  error  for  the  lower  noise  values  (^<0.005 
radian)  and  deviate  significantly  for  anr 0*018 
radian.  (See  Figures  4  and  5.) 

c.  The  variation  of  RMS  residual  with  measurement 
noise  is  essentially  linear  with  a  slope  of  unity 
(Figure  6) .  This  indicates  that  the  RMS  residual, 
as  calculated  by  the  coefficient  extraction  pro¬ 
gram,  is  a  good  estimate  of  the  amount  of  measure¬ 
ment  noise  in  the  data  when  this  is  the  only 

type  of  noise  present. 

d.  All  variations  of  the  normalized  estimated  standard 
deviations  of  the  extracted  coefficients,  crjn, 
increase  linearly  with  measurement  error  (Figures 

7  through  10) .  Approximate  values  of  the  slopes 
of  these  variations  are  given  in  Table  V.  These 
normalized  uncertainty  ratios  give  the  relative  (to 
each  other)  uncertainty  which  can  be  expected  when 
extracting  coefficients  from  data  containing 
measurement  noise  using  the  given  coefficient 
extraction  program. 

e.  The  true  value  of  a  given  coefficient  is  contained 
within  the  interval  defined  by  its  extracted 
value  ±3aj  for  every  coefficient  extracted  from 
the  nine  sets  of  data. 
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For  this  portion  of  the  analysis,  attempts  were  made 
to  extract  aerodynamic  coefficients  from  twelve  basic 
sets  of  dynamic  data.  Each  data  set  was  made  up  of  2G1 
discrete  points  resulting  from  integrating  the  equation 
of  motion  in  increments  of  0.005  second  for  4.00  seconds; 
the  pitch  angle  was  output  every  four  integration  steps 
so  that  the  time  between  data  points  was  0.020  second. 

The  time  between  data  points  used  in  this  portion  of  the 
analysis  is  different  from  that  previously  used.  This 
change  was  made  so  that  the  time  between  data  poix  ts  would 
correspond  to  that  used  by  Turner^  j.n  actual  experiment¬ 
ation. 

Of  the  twelve  basic  data  sets,  nine  contain  system 
noise  and  measurement  errors ,  whereas  three  contain  only 
measurement  errors.  Coefficients  were  extracted  from  these 
latter  three  to  determine  if  the  above-mentioned  change  in 
the  time  increment  between  data  points  had  any  appreciable 
effect  on  the  extracted  coefficients  or  their  uncertainties. 
The  estimates  of  the  uncertainties  were  improved  due  pri¬ 
marily  to  the  fact  that  a  longer  data  record  was  being 
analyzed. 

The  desired  mean  values  of  the  frees tream  velocity 
vector  direction  fluctuations  and  measurement  er  -ors  for 
all  data  sets  were  zero.  The  desired  standard  de  'iation 
of  the  measurement  errors,  ©m,  was  0.00582  radian 
(3crm  =  1.0  degree)  for  all  sets;  the  desired  standard 
deviation  of  the  velocity  magnitude  fluctuations,  ov, 
was  5.0  ft/sec  for  ail  nine  data  sets  containing  system 
noise.  The  desired  standard  deviations  of  the  velocity 
direction,  <ra,  were  0.02910  radian  (3aa  =  5  degrees), 

0.05820  radian  (30a  =  10  degrees),  and  0.11640  radian 
(3 oa  =  20  degrees) ,  with  three  different  sets  of  data 
being  generated  at  each  noise  level. 

The  random  velocity  fluctuations  act  as  a  forcing 
function  for  the  equation  of  motion  and  cause  oscillations 
even  if  the  vehicle  has  no  initial  displacement.  The 
resultant  maximum  amplitudes  of  these  forced  oscillations 
are  usually  within  a  certain  magnitude  or  noise  band  width. 
The  widths  of  the  noise  bands  for  the  data  used  in  this 
analysis  are  generally  equivalent  to  the  30 a  values,  and 
these  have  been  used  to  calculate  the  percent  system  noise 
values.  The  approximate  average  peak  amplitude  used  in  the 
percent  noise  calculations  is  the  mean  positive  peak  ampli¬ 
tude  for  the  10  cycles  that  result  when  the  nominal  equation 
of  motion  is  integrated  for  approximately  four  seconds. 


SUNMARY  OF  RESULTS  FOR  EXTRACTING  AERODYNAMIC  COEFFICIENTS  FROM  DYNAMIC  DATA  CONTAINING 


Figure  11  is  an  example  of  one  of  the  program  fits 
to  a  noisy  data  set.  Figures  12  through  15  show  the  varia¬ 
tion  of  percent  error  for  the  extracted  coefficients  with 
system  noise.  The  variations  of  the  normalized  estimated 
standard  deviations  of  the  extracted  coefficients  with 
system  noise  level  are  shown  on  Figures  16  through  19. 

An  analysis  of  Table  VI  and  Figures  12  through  19 
reveals  the  following: 

a.  The  errors  in  the  linear  static  pitching  moment 
coefficient  derivative  parameter,  C3,  are  gen¬ 
erally  less  than  5  percent  for  the  entire  range 
of  system  noise  considered.  The  errors  in  the 
corresponding  nonlinear  term,  C5,  are  generally 
less  than  10  percent.  (See  Figures  12  and  13.) 

b.  The  pitch  damping  parameters  show  significant 
error  for  all  nonzero  system  noise  levels.  The 
extracted  values  of  the  nonlinear  term,  Cg ,  are 
sometimes  in  error  by  several  orders  of  magnitude 
at  the  higher  noise  levels.  (See  Figures  14 

and  15.) 

c.  The  estimated  standard  deviations  of  the  extracted 
coefficients  generally  increase  as  the  system  noise 
increases.  The ^estimated  standard  deviations 

of  coefficients  "extracted  from  dynamic  data 
containing  both  system  noise  and  measurement 
errors  are  generally  too  small  and  do  not  reflect 
the  true  uncertainty  of  the  extracted  coeffi¬ 
cients  as  was  true  in  the  previous  case  when 
the  data  contained  only  measurement  errors. 

(See  Table  VI  and  Figures  16  through  19.) 

d.  All  attempts  at  extracting  the  complete  set  of 
coefficients  from  data  with  a  system  noise  band 
of  approximately  20  degrees  failed  (see  Table  VI) . 
The  failures  resulted  when  a  value  of  Cg  was 
eventually  calculated  by  the  iterative  process 
which  caused  the  solution  to  the  equation  of 
motion  to  diverge. 

e.  The  divergence  problem  can  sometimes  be  circum¬ 
vented  by  not  attempting  to  extract  Cmq2  prom 
extremely  noisy  data.  The  resulting  coefficients 
that  are  extracted,  Cmao/  Cm'-*2/  an$  Ongo/  have 
accuracies  comparable  to  those  extracted  from 
noisy  data  with  a  ±10  degrees  system  noise  band. 
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Program  Fit:  C  =  -1.97  -  24.7a 
ma 


Figure  1.  Example  of  Program  Fit  to  Data  Containing  Measurement  Errors 
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Deviation  of  C6  with  System  Noise 


SECTION  IV 


DEVELOPMENT  OF  THE  EXTENDED  KALMAN  FILTER  FOR 
ESTIMATING  PARAMETERS  AND  THEIR  UNCERTAINTIES 


1.  Introduction 


This  chapter  presents  an  abridged  derivation  of  the 
Kalman  filter  for  discrete  and  continuous  linear  systems, 
followed  by  a  general  statement  of  the  extended  Kalman 
filter  for  continuous  nonlinear  systems.  The  extended 
filter  is  then  used  as  a  base  for  the  development  of  a 
specific  algorithm  for  estimating  states  and  parameters 
of  the  second  order  equation  of  pitching  motion  for  a 
missile  being  forced  by  random  disturbances. 


2.  Development  of  the  Kalman  Filter 

The  original  derivation  o£  the  Kalman  filter  was  presented  by 
Kalman18  in  1960 -for  multistage  systems  making  discrete  linear  transi¬ 
tions  from  one  stage  to  another.  Kalman  and  Bucy19  gave  the  analogous 
development  for  continuous  linear  systems  approximately  one  year  after 
the  first  work  was  published.  The  purpose  of  the  linear  filter  is  to 
provide  estimates  of  the  state  of  a  system  by  making  use  of  measure¬ 
ments  of  all,  or  some,  of  the  state  vector  components  of  the  system. 

The  system  is  assumed  to  be  operating  in  the  presence  of  random  distur¬ 
bances,  the  statistical  properties  (i.e.,  mean  and  variance)  of  which 
are  known.  The  measurements  of  the  state  vector  components  of  the 
system  are  also  assumed  to  have  randan  errors  of  known  statistics. 

In  addition  to  the  original  derivations  of  Kalman  and  Bucy, 
several  other  methods  offering  various  degrees  of  insight  but  leading 
to  the  same  results  are  available.  A  brief  development  taken  primarily 
from  Bryson  and  Ho20  is  presented  here.  Other  derivations  or  develop¬ 
ments  of  the  filter  equations  are  given  by  Jazwinsky,3 1  Sorenson,  2 
and  Barham.3 3 

The  treatment  given  here  starts  with  a  static  system  and  is 
extended  for  a  single-stage  linear  transition  which  leads  directly  to 
linear  multistage  process.  Finally,  by  making  use  of  a  limiting 
process,  the  desired  form  of  the  equations  for  a  continuous  linear 
dynamic  system  is  given. 

a.  Static  System 

The  problem  at  hand  is  to  estimate  the  n-component  state  vector 
X  of  a  static  system  using  the  p-component  measurement  vector,  z,  con¬ 
taining  random  errors,  v,  which  are  independent  of  the  state.  The 
measurement  vector  can  be  represented  as 
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z  ---  HX  +  V 


i 


(21) 


where  H  is  a  known  p  x  n  matrix.  Conditions  on  the  measure¬ 
ment  errors  are 


E  (v)  =  0 


(22) 


E(WT)  -  R 


(23) 


where  R  is  a  known  matrix  of  dimension  p.  It  is  assumed 
that  a  prior  estimate  of  the  state,  designated  as  X,  is 
available  and  also  that  the  covariance  of  the  prior  esti¬ 
mate,  M,  is  known.  Thus 


M  =  E  [  (X-X)  (X-X)  T] 


(24) 


where  M  is  of  dimension  n. 

The  desired  estimate  of  X,  taking  into  account  the 
measurement,  z,  is  the  weighted-least-squares  estimate, 
ft,  which  minimizes 


J  =  t(X-X)T  M-1  (X-X)  +  (z-HX)T  R"1  (z-HX)].  (25) 


Differentiating  Equation  (2.5)  with  respect  to  X,  setting 
the  resulting  Equation  to  zero,  and  solving  for  X  yields 


ft  -  X.  +  PHTR”‘*  (z-HX) 


(26) 


where 


P-1  -  HF1  +  HTR~1H 


(27) 
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The  quantity  P  is  the  covariance  matrix  of  the  error  in 
the  state  estimate  after  measurement,  X,  i.e,. 


P  =  E[  (X-X)  (X-X)T]  .  (28) 
b.  Single-Stage  Transitions 

It  is  now  desired  to  estimate  the  state  of  a  system 
that  makes  a  discrete  transition  from  state  0  to  state  1 
according  to  the  linear  relation 


Xi  =  4>X  +  JT  w 

1  o  o  o  o 


(29) 


where  is  a  known  transition  matrix  of  dimension  n  and 
rQ  is  a°known  n  x  r  matrix.  The  mean  and  variance  of  the 
random  forcing  vector  are  assumed  to  be  known  and  are 
given  by 


E(wo)  = 


w 


El(w0-S0)  (W0-W0)T]  =  Qc 


(30) 


The  statistical  properties  of  the  random  state  vector  are 
assumed  to  be  known  initially  as 


®<V  '  x0  '  E[<VV  <Vxo>Tl  =  po  •  (31) 


It  is  also  assumed  that  XQ  and  wQ  are  independent.  Prom 
this  inrormation  Xi  is  a  random  Sector  whose  mean  value 
and  covariance  are 


X  = 


M 


1  - 


6  x  +  r  w 
o  o  o  o 

(32) 

T  T 

6  p  4,  +  r  q  r 
yO  OyO  0*0  o 

(33) 

Suppose  now  that  measurements  of  the  state  are  made  after 
the  transition  to  state  1;  then  from  Equations  (26)  and 
(27)  the  best  estimate  of  Xj  is 

/t i  =  Xl  -r  P^fR-5  (z-HjXj )  ,  (34) 
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where 


Px  =  (Mf1  +  H^r-Ihj)  ~l 

=  Mx  -  MxHi  +  P.J1)  HjMx 


(35) 


c.  Linear  Multistage  Processes 

The  developments  of  the  two  preceding  subsections  can 
be  extended  for  linear,  stochastic,  multistage  processes. 
Given  the  following  difference  equation  system  model 


=  4>-X.  +  r.w. 

i+l  ii  ii 


(36) 


where 

E<V  =  ?o 

E  (Wi)  =  Wi 

E[(X0-X0)(X0-X0)T1  =mo 
E [  (w^-wi)  (w^j)]  -  QjSy 

E  t (X0-X0)T]  -  0 
and  measurements  of  the  state 


(37) 

(38) 

(39) 

(40) 

(41) 

(42) 
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Equations  (46) ,  (47)  and  (48)  are  the  discrete  Kalman 
filter  with  the  state  variable  covariances  given  by 
Equations  {'’9)  and  (50).  As  can  be  seen  from  the  above 
equations,  the  Kalman  filter  is  essentially  the  same  as 
the  system  model  [Equation  (36)).  The  differences  are 
(1)  the  actual  system  noise,  which  is  random  from  one  stage 
to  the  next  in  Equation  (36) ,  is  replaced  by  its  mean  or 
expected  value,  and  (2)  there  is  a  correction  term  based 
on  the  difference  between  the  actual  measurement  of  the 
state  and  its  predicted  value i  The  difference  term  is 
multiplied  by  a  gain,  K^,  which  is  essentially  the  ratio 
of  the  uncertainty  in  the  state  to  the  uncertainty  in 
the  measurement.  If  the  covariances  of  the  measurement 
errors  are  large,  the  gain  will  be  relatively  small  and 
the  corresponding  difference  term  will  have  little  effect 
on  the  estimate  of  the  state;  if,  however,  the  system 
noise  is  relatively  large  or  the  measurement  errors  are 
very  small,  the  gain  will  be  large  and  differences  between 
the  actual  and  predicted  measurements  of  the  state  at  a 
given  stage  have  increased  significance  in  the  state 
variable  estimates. 

d.  The  Continuous  Kalman  Filter 

By  applying  a  limiting  process  to  the  discrete  filter 
with  the  time  between  stages  tending  to  zero,  the  linear 
system  model  becomes 


4 

X  =  F(t)X  +  G  (t)w(t) 


(51) 


and  the  continuous  Kalman  filter  is  given  by 


X  *»  FX  +  Gw  +  PHTR_1  (z-HX)  (52) 

P  -  FP  4  PFT  +  GQGT  ~  PK^R“ 1 HP  .  (53) 


3.  The  Extended  Kalman  Filter 

The  extension  of  the  Kalman  filter  for  estimating 
the  states  of  nonlinear  systems  in  the  presence  of  noise 
has  been  given  by  Erysoh  and  Ho^O  and  Jazwinski ,31  among 
others ,  as 


41 


X  =  f  (X,  t)  +  G(t)w(t)  +  p[§]  R"1 

[ z (t) -h  (X,t)  ] 

(54) 

P  =  —  P  +  p  1  —1  ~  +  gqgt  -  P  T— 1 
ax  Laxj  bxj 

R-i  rihi  P 

L  3XJ  e 

(55) 

for  the  nonlinear  system 

X  =  f(X,t)  +  G(t)w(t) 

(56) 

with  measurements 

z(t)  =  h(X,t)  +  v(t) 

(57) 

where 

E[w(t)  ]  =  w  (t) 

(58) 

E  [w  ( t)  -w  (t)  ]  [w  (t 1 )  ]  =  Q  { t)  6  (t-t ' ) 

(59) 

E[x(ton  =  xo 

(60) 

E(X(t0)-X0][X(t0)-X0]T  =  PQ 

(61) 

EiX(t0)-X0’J  [w(t)-w(t)]T  =  0 

(62) 

E[v(t)]  =  0 

(63) 

E  [v  (t)  vT  (t ' )  ]  =  R  (t)  6  (t~t 1 ) 

4 

(64) 

42 


0 


(65) 


E[X(to)-XQ]  [v(t)]T  = 


E[w(t)-w(t)  ]  [v(t')]T  =  0 


(66) 


4.  The  Parameter  Estimation  Algorithm 

In  this  subsection  a  specific  parameter  and  state 
estimation  algorithm  using  the  extended  Kalman  filter  is 
developed  for  the  basic  equation  of  pitching  motion  pre¬ 
viously  analyzed  with  the  Chapman-Kirk  technique. 

The  nonlinear  equation  of  motion,  or  system  model, 
is  thus 


d  +  (C4+Cgci2)a  +  (C3+C502)a  =  w(t) 


(67) 


where 


E  [w  (t)  ]  =  0 


E[w(t)w(t')  ]  =  q6(t-t')  .  (68) 


The  measurements  of  the  state  of  the  system  are  assumed 
to  be  given  by 


z  (t)  =  a  (t)  +  v  (t) 


(69) 


where 


E  (v  (t)  ]  =  0 


E[v(t)vT(t')  ]  =  r6  (t-t') 


(70) 


Now  to  reduce  Equation  (.67)  to  the  required  system  of 
first  order  differential  equations,  let 
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ct  =  X ! 


a  =  X2  .  (71) 

At  this  point,  to  estimate  the  aerodynamic  parameters 
appearing  in  the  equation  of  motion  in  addition  to  the 
state  variables ,  the  state  vector  is  augmented  by  setting 

C3  =  X3 
C4  =  X4 
C5  =  X5 

C6  =  X6  (72) 


with  the  constraints 


♦ 

x3  =  0 
x4  =  0 

X5  =  0 

X6  =  0  .  (73) 


Making  use  of  Equations  (71)  ,  (72)  ,  and  (73)  ,  Equation  (67) 
now  reduces  to  the  following  nonlinear  system  of  first 
order  equations : 
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X2  =  -(X4+X6X12)X2-(X3+X5Xi2)X1  +  w(t) 

X3  =  0 
X4  =  0 


AS  =  o 

x6  =  0  ,  (74) 


with  linear  measurements 


z  =  Xi  +  v 


(75) 


Comparing  Equations  (74)  and  (75)  with  Equations  (56)  and 
(57) ,  the  following  matrices  may  be  identified: 


G  ( t)  = 


0 

H  (t)  = 

Xl 

1 

0 

0 

0 

0 

0 

0 

0 

0 

y 

0 

f(X,t)  = 


X2 

-(X4+X6X12)X2  -  (X3+X5X!2)X3 

0 

0 

0 

0 


.  (76) 
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Now,  since  the  mean  value  of  the  system  noise  is  assumed 
to  be  zero  [Equation  (68)]  and  because  the  measurements 
of  the  state  are  related  linearly  to  the  state  [Equation 
(75)],  the  extended  Kalman  filter,  previously  given  by 
Equations  (54)  and  (55)  ,  can  be  simplified  somewhat  to 


x  =  (X,t)  +  PHT  \l\  [z-HX] 

L-J 


P  =  —  P 

3X 


+  P 


3f 

T 

„  „T 

fll 

ax 

+  GqG  -  PH 

HP 


where 


H  =  [100000] 


Applying  Equations  (76)  and  (77)  yields 


*1 

A 

x2 

X3 

A 

x4 

X5 

x6 


(Xli+X6Xl2)X2  -  (X3+X5X12)X1 
0 
0 
0 
0 


P11 

P  1  2 

P 1  3 

P 1  4 

P 1  5 

pl  6 

P  1  2 

P2  2 

p2  3 

p2  4 

P2  5 

p2  6 

P  1  3 

P2  3 

P  3  3 

P  3  4 

P  3  5 

P  3  6 

P14 

P  2  4 

p  3  4 

P4  4 

P4  5 

P4  6 

p  1  5 

p  2  5 

P  3  5 

P4  5 

P  5  5 

P  5  6 

±(z-[l  00000] 


Pi  6  P26  p36  P46  P 5 6  p66 


0 


Carrying  out  the  prescribed  multiplications  and  equating 
like  components  of  the  resulting  two  column  matrices 
yields  the  desired  filter  equations, 


Xi  =  x2  +  ^  (Z-Xj) 

X2  =  -(X4+X6X^)X2  -  (Xa+XgX^Xj  +  -LZ.  (z-Xj) 

:  p 

X3  =  (Z-Xl ) 

:  p 

X4  -  "•  (z-Xj) 

;  P 

X5  =  (z-Xi) 


(79) 


The  necessary  covariance  equations  are  obtained  from  the 
matrix  Ricatti  equation  [Equation  (78)]  making  use  of 
Equations  (76)  and  the  fact  that 


0 

-2X1X2X6-X3-3XfX5 

0 

0 

0 

0 


1 

~(X4+X6X2) 

0 

0 

0 

0 


47 


(80) 


The  results  of  substituting  Equations  (76)  and  (80)  into 
(78) ,  carrying  out  the  rather  tedious  matrix  multipli¬ 
cations  ,  and  equating  like  terms  are 

PU  =  2P12-|pn2 

P 1  2  =  P22“aPii-BP12-XiPi  3“X2Pi4'*CPi5-DPi6-pPiiPi2 
P 1 3  =  P23“^PllP13 
P 1  4  =  P  2  4  “  jP  1 1 P 1  4 
P 1  5  =  P  2  5  ~*^P  1 1 P 1  5 

p  1  6  =  P  2  6  — pP 1  1 P 1  6 

P2  2  =  q-2 (AP12+EP22+X1P23+X2P24+CP2  5+DP26) "jPl22 

P23  =  -AP13-BP23-X1P33-X2P34'-CP35-DP36“Pi2Pi3 
P2i,  =  — AP i  4— BP 2 4**XlP 3 4~X2P 4 4 *"CP 4  5— PPh g*~P  1  2P 1  4 
P25  =  -AP 1 5 ~BP 2  5-X1P 3  5~X2P 4  5 -CP 5 3 “DP 5 5— jP 1 2P 1 5 
P26  =  “AP 3  g -BP 2 e -X 1 P 3  6 *~X2P 4  6 -<-P 5  5 -DP 5 g “pP  1  2P 1  6 
P33  =  "^Pl32 
P34  =  "P13P14 

P35  *  -7P13P15 

P  3  6  =  "Pl3Pl6 
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P44  =  “^Pl42 
P45  =  4P14P15 
P46  “  “£Pl4Pl6 
P55  =  ~^P 1 5  2 

*  t 

PS  6  =  “jPl5Pl6 
P66  =  -jPl62 

where 


A  =  2X1X2^6  +  X3  +  3X^X5 
B  =  X4  +  X6Xf 

C  =  X^ 

D  =  X*X?  .  (82) 


The  desired  estimates  of  the  states  and  parameters 
are  obtained  by  numerically  integrating  the  filter  equations 
[Equations  (79)]  and  their  covariances  [Equations  (81)]. 
Initial  estimates  of  the  states  and  their  covariances  as 
well  as  the  variances  of  the  system  noise  and  measurement 
errors  are  assumed  to  be  available.  In  the  event  that  the 
data  consist  of  discrete  measurements,  which  is  usually 
the  case,  che  constant  time  between  data  points  should  be 
equa.l  to  the  numerical  integration  step  size.  This  is 
necessary  because  the  state  estimates  are  updated  at  each 
integration  increment  and,  in  doing  this,  the  filter  requires 
knowledge  of  the  difference  between  measured  and  estimated 
state  values. 

If  the  filter  successfully  adjusts  the  true  states  and 
the  imbedded  parameters  so  that  Xj  is  an  adequate  match  to 
the  data,  z,  the  time  derivatives  of  the  filter  equations 
for  the  parameters  become  small  [see  Equations  (79)]  and 
the  parameters  reach  steady-state  or  near  steady-state 
conditions . 
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SECTION  V 


USE  OF  THE  EXTENDED  KALMAN  FILTER  FOR 
ESTIMATING  PARAMETERS  AND  THEIR  UNCERTAINTIES 


1.  Introduction 


This  section  presents  an  analysis  using  the  extended 
Kalman  filter  with  parameter  augmented  state  vector  to 
determine  aerodynamic  coefficients  and  their  uncertainties 
from  noisy  dynamic  data.  A  description  of  the  digital 
computer  program  used  in  the  analysis  is  given  first, 
followed  by  results  obtained  with  this  program  from  a 
variety  of  noisy  data  sets. 

2 .  The  Extended  Kalman  Filter  Program 

The  extended  Kalman  filter  program  is  written  in 
Fortran  IV  for  use  primarily  on  an  IBM  model  360/65 
digital  computer.  The  basic  function  of  the  program  is 
to  integrate  numerically  the  27  first  order  differential 
equations  which  are  given  in  the  previous  chapter  and 
which  comprise  the  extended  Kalman  filter  with  parameter 
augmented  state  vector.  The  specific  functions  of  the 
main  program  and  its  subroutines,  as  well  as  the  required 
input  data,  are  described  in  the  following  paragraphs. 

A  listing  of  the  complete  program  is  given  in  Appendix  II. 

a.  Main  Program  Functions 

The  main  program  reads  and  writes  the  input  data 
and  calls  the  numerical  integration  subroutine, 
ADDUM. 

b .  Subroutine  Functions 

The  names  and  functions  of  each  of  the  three  sub¬ 
routines  that  are  used  in  the  extended  Kalman  filter 
program  are  given  below. 

(1)  ADDUM 

Subroutine  ADDUM  integrates  the  27  filter  and 
covariance  equations  using  a  fourth  order  Runge- 
Kutta  method  for  starting,  followed  by  a  fourth 
order  Adams-Bashforth  predictor-corrector  method 
for  running.  It  is  essentially  identical  to  the 
subroutine  of  the  same  name  used  in  the  Chapmans 
Kirk  coefficient  extraction  program  and  is  described 
in  detail  in  Reference  26 .  Subroutines  XDOT  and  OUT 
are  called  from  ADDUM. 
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(2) 


XDOT 


Subroutine  XDOT  computes  values  of  the  time  deri¬ 
vatives  for  the  differential  equations  being 
integrated  by  ADDUM. 

(3)  OUT 

This  subroutine  writes  the  output  of  the  numerical 
integration . 

*- *  Hsqurred  Input  Data 

The  program  reads  four  categories  of  input  data. 
These  categories  and  the  specific  data  in  each  are 
delineated  in  Appendix  IX .  The  format  and  units 
of  the  entries  on  a  specific  data  card  can  be 
determined  from  the  program  listing  and  nomenclature 
list,  which  are  also  given  in  Appendix  II. 

3 .  Analysis  of  the  Extended  Kalman  Filter 

The  use  of  the  extended  Kalman  filter  for  estimating 
the  parameters  of  interest  is  analyzed  for  linear  and  non¬ 
linear  systems  with  data  containing  only  measurement  errors 
as  well  as  data  containing  both  system  noise  and  measurement 
errors.  The  analysis  begins  with  a  linear  system  and  data 
containing  only  measurement  errors  and  progresses  through 
increasing  stages  of  difficulty  up  to  nonlinear  systems  and 
data  which  have  been  corrupted  by  both  system  noise  and 
measurement  errors.  The  previously  discussed  computer 
program  and  algorithm  require  no  modifications  to  consider 
the  linear  case;  by  initially  setting  the  parameters  that 
are  the  numerical  coefficients  of  the  nonlinear  terms  in 
the  equation  of  motion  (and  their  variances)  equal  to  zero, 
the  extended  filter  for  a  nonlinear  system  reduces  to  the 
one  required  for  a  linear  system. 

All  data  used  in  the  analysis  were  generated  by  the 
previously  mentioned  computer  program,  UFNOISE  (Reference 
29) .  The  true  values  of  the  constants  used  in  the  data 
generation  are  those  previously  given  in  Table  III  with  the 
exception  of  the  initial  value  of  the  pitch  angle  for  the 
linear  system.  This  initial  displacement  is  a  more  realistic 
0.1745  radian  for  the  linear  cases  considered. 

a.  Linear  Systems  with  Measurement  Errors  Only 

After  the  extended  Kalman  filter  computer  program 
had  been  constructed  and  checked,  it  became  apparent 
that  the  technique  would  probably  be  best  under¬ 
stood  by  considering  relatively  simple  cases  at  first 
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and  then  progressing  to  more  difficult  ones  as 
confidence  in  the  technique  was  gained.  To  this 
end,  the  first  success  with  the  program  was  rea¬ 
lized  when  analyzing  data  containing  only  measurement 
errors  for  a  linear  system. 

The  pitch  angle  data  used  in  this  part  of  the 
analysis  were  generated  bv  integrating  the  equation 
of  motion 


C40  +  C  3  a  =  0 


(83) 


and  superimposing  random  Gaussian  errors  on  the 
output.  The  standard  deviation  of  the  errors  is 
0.00588  radian,  which  corresponds  to  3 a  measurement 
errors  of  approximately  1.0  degree.  The  numerical 
integration  step  size  used  in  generating  the  data 
was  0.005  second,  and  the  pitch  angle  was  output 
eveiy  integration  increment.  The  equation  was  in¬ 
tegrated  for  a  total  of  2.00  seconds. 

(1)  Basic  Filter  Performance 

The  results  of  using  the  extended  Kalman  filter  to 
identify  the  correct  values  of  the  parameters  of 
interest,  X3  and  X4  ,  are  shown  on  Figures  20  and 
21.  These  figures  show  time  histories  of  the  per¬ 
cent  errors  in  the  estimated  values  of  the  parameters 
that  were  computed  by  the  extended  filter  program 
when  analyzing  the  noisy  data  described  previously. 
The  percent  error  is  defined  by 


(X.-X.)'IOO 

PERCENT  ERROR  (X.)  =  - L-rH - 

3  xj 


(84) 


The  input  for  the  measurement  error  variance  was  the 
one  previously  computed  by  the  simulation  program 
based  on  the  errors  that  were  actually  put  into  the 
data.  The  initial  values  of  the  parameter  variances 
were  computed  from  knowledge  of  the  true  value  of 
the  parameters  and  the  arbitrarily  selected  initial 
values  of  the  parameters.  These  sample  variances 
are  defined  by 


p..(0)  =  [X_.  (0)  -X .  (0)  ]  2 
33  J  3 


(85) 
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This  method  of  initializing  the  parameter  variances 
was  chosen  merely  to  generate  the  required  initial 
numerical  values.  In  actual  experimentation,  the 
initial  variance  values  would  depend  on  the  method 
of  selecting  the  initial  parameter  estimates  and 
prior  knowledge  as  to  hov  accurately  these  initial 
parameter  estimates  were  in  relation  to  the  true 
values  of  the  parameters. 

The  initial  values  of  all  covariances  were  chosen 
as  zero.  This  implies  that  errors  in  the  estimates 
of  the  individual  state  vector  components  are  initially 
uncoxreletted.  The  initial  values  of  the  variances 
for  the  pitch  angle  and  pitch  angle  rate  were  also 
chosen  as  zero.  For  wind  tunnel  dynamic  experimen¬ 
tation  where  the  model  is  initially  held  rigid  at 
some  given  displacement  to  the  wind  stream,  this 
seems  to  be  a  valid  assumption. 

As  can  be  seen  from  Figures  20  and  21,  the  filter 
does  an  excellent  job  of  identifying  the  parameters 
of  interest.  The  X3  parameter,  which  corresponds 
to  the  Cmct  term  and  which  is  initially  in  error  by 
25  percent,  is  identified  with  approximately  zero 
percent  error  in  less  than  0.3  second.  The  correct 
identification  of  the  damping  parameter  takes  slightly 
more  time  but  the  results  are  of  equal  quality. 

(2)  Effects  of  Variations  in  the  Initial  Parameter  Variances 

This  subsection  presents  some  effects  on  the  near 
steady-state  estimates  of  the  parameters  of  interest 
and  rheir  near  steady-state  variances  for  various 
values  of  the  initial  parameter  variances.  The 
parameter  variances  were  initialized  to  values  that 
were  25  percent  higher  than  those  used  in  the  pre¬ 
vious  sub  section  and  also  to  values  that  were  25 
percent  lower  than  the  referenced  values.  The  effects 
of  these  initializations  are  shown  in  Figures  22 
through  25  and  in  Table  VII  for  the  linear  system 
of  interest. 

Table  VII  is  a  summary  of  the  percent  errors  in  the 
parameter  variances  and  the  normalized  parameter 
uncertainties  for  the  above-mentioned  variations  in 
the  initial  parameter  variances. 

The  same  results  are  depicted  graphically  in  Figures 
22  and  23.  The  numerical  results  in  Table  VII  are 
based  on  the  near  steadyr-state  parameter  and  variance 
values  that  result  after  the  filter  has  integrated 
for  1.5  seconds.  These  correspond  to  the  last  point 
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TABLE  VII.  EFFECTS  OF  VARIATIONS  IN  THE  INITIAL 
PARAMETER  VARIANCES  ON  NEAR  STEADY- 
STATE  PARAMETERS  AND  THEIR  VARIANCES 
(LINEAR  SYSTEM,  MEASUREMENT  ERRORS  ONLY) 


Variation  in 
Initial  Parameter 
Variance* 
(Percent) 

Errors  in  Near 
Steady-State 
Parameters 
(Percent) 

Normalized  Near 
Steady-State 
Variances 
(Percent) 

A 

X3 

A 

x4 

*^3  3 

n 

/P44 

n 

-25 

0.17 

0.55 

0.12 

0.91 

0 

C.I9 

0.46 

0.11 

0.91 

25 

0.22 

0.38 

0.11 

0.91 

♦relative  to  sample  variance 


shown  on  the  supporting  figures.  The  normalized 
parameter  uncertainties  are  defined  by 


/PTT 

PERCENT  (v'FTT)  =  — 11  .  100  .  (86) 

33  n 


Time  histories  of  the  parameter  variances,  P33  and 
P4 4 ,  are  shown  in  Figures  24  and  25.  All  of  the 
information  presented  in  these  figures  indicates 
that  initial  values  of  parameter  variances  in  the 
range  investigated  have  no  effect  on  final  parameter 
estimates  or  their  variances  when  using  the  extended 
Kalman  filter  to  analyze  data  for  a  linear  system 
containing  only  measurement  errors. 

b .  Linear  Systems  with  Measurement  Errors  and  System 
Noise 


Data  for  this  portion  of  the  analysis  were  generated 
with  the  computer  program  UFNOISE  by  numerically 
integrating  the  following  equation  of  motion 


a  +  C40  +  C3C1  =  w(t) 


(87) 


and  superimposing  random  Gaussian  errors  on  the  out¬ 
put.  The  forcing  function,  w(t) ,  is  also  random 
in  nature  and  Gaussian  with  zero  mean.  The  standard 
deviation  of  the  system  noise  was  4.80  rad/sec2, 
which  is  of  sufficient  magnitude  to  drive  the 
oscillations  in  a  noise  band  of  approximately  0.0872 
radian  (5  degrees)  for  zero  initial  displacement. 

The  standard  deviation  of  the  measurement  errors 
was  0.00556  radian.  Also,  as  before,  the  equation 
of  motion  was  integrated  in  increments  of  0.005 
second  for  a  total  of  2.00  seconds. 

(1)  Basic  Filter  Performance 

The  basic  performance  of  the  filter  is  demonstrated 
in  Figures  .  •>  and  27  which  show  time  histories  of 
the  percent  error  in  X3  and  X4 .  As  can  be  seen  from 
these  figures,  the  initial  estimates  of  X3  and  X4 
are  25  percent  in  error  and  these  estimates  are 


corrected  to  within  approximately  1  percent  and 
3  percent  of  their  respective  true  values.  The 
noise  variances  used  are  those  given  by  the 
simulation  which  generated  the  pitch  angle  data. 

The  initial  parameter  variances  are  calculated 
from  Equation  (85) . 

(2)  Effects  of  Variations  in  the  Initial  Parameter 
Variances 

The  sensitivity  of  the  filter  to  variations  in 
the  initial  parameter  variances  is  given  in 
Table  VIII  and  in  Figures  28  through  32.  Table 
VIII  is  a  summary  of  the  percent  error  in  the 
parameters  of  interest  and  their  normalized  un¬ 
certainties  for  variations  in  the  initial  para¬ 
meter  variances  of  -25  percent,  0.0  percent  and 
25  percent  relative  to  the  sample  variances. 

Figures  28  and  29  show  the  variation  of  the  near 
steady-state  errors  in  the  parameters  and  their 
normalized  uncertainties  as  functions  of  the  per¬ 
cent  variations  in  the  initial  parameter  variances. 
Figure  30  shows  the  two  time  histories  of  the 
error  in  the  estimate  of  X4  that  result  for  two 
different  initial  values  of  the  parameter  variances. 
One  of  the  initial  variances  is  the  sample  variance 
based  on  the  initial  and  true  values  of  the  para¬ 
meter  (Equation  (85) ]  ,  and  the  other  is  25  percent 
less  than  this  value. 

The  X3  error  time  histories  are  essentially  iden¬ 
tical  for  both  cases  and  are  not  presented.  The 
variations  of  the  parameter  variances  with  time, 
for  the  three  initial  values  considered,  are  shown 
in  Figures  31  and  32. 

From  the  information  given  in  the  figures  mentioned 
above,  it  is  obvious  that  neither  X3  nor  P33  are 
affected  by  variations  in  the  initial  variances 
within  the  ±25  percent  range  considered.  The  es¬ 
timates  of  the  damping  parameter  as  well  as  its 
variance  are  affected,  however,  by  these  variations. 
The  near  steady-state  damping  pararnt  :.er  estimates 
vary  almost  linearly  with  initial  variance  values 
from  a  low  of  0.47  percent  to  a  high  of  4.95  per¬ 
cent  (see  Figure  28) ,  The  variation  in  the  normalized 
uncertainty  is  less  pronounced,  ranging  from  a  low 
of  14.39  percent  to  a  high  of  15.81  percent  (see 
Figure  39).  Nevertheless,  the  near  steady-state 
uncertainties  are  of  sufficient  magnitude  so  that 
the  true  values  of  the  coefficients  are  within  less 
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TABLE  VIII.  EFFECTS  OF  VARIATIONS  IN  THE  INITIAL 
PARAMETER  VARIANCES  ON  NEAR  STEADY- 
STATE  PARAMETERS  AND  THEIR  VARIANCES 
(LINEAR  SYSTEM,  MEASUREMENT  ERRORS 
AND  SYSTEM  NOISE) 


Variation  in 
Initial  Parameter 
Variance* 
(Percent) 

Errors  in  Near 
Steady- State 
Parameters 
(Percent) 

Normalized  Near 
Steady-State 
Variances 
(Percent) 

/s 

x3 

A 

•^33 

n 

n 

-25 

0.72 

-0.47 

2.42 

14.39 

0 

0.38 

2.66 

2.42 

15.18 

25 

0.98 

4.95 

2.42 

15.81 

♦relative  to  sample  variance 

than  one  standard  deviation  of  their  estimated 
value  regardless  of  the  initial  parameter  variance, 
when  the  standard  deviation  is  taken  as  the  square 
root  of  the  near  steady-state  parameter  variance. 
Attention  is  also  called  to  the  fact  that  since 
system  noise  is  now  present  in  the  data,  the  near 
steady-state  parameter  variances  are  of  greater 
magnitude  than  was  the  case  when  only  measurement 
errors  were  present  (see  Figures  24,  25,  31,  and 
32)  . 

(3)  Filter  Response  to  Large  Errors  in  the  Initial 
Parameter  Estimates 

The  results  of  using  the  extended  Kalman  filter  to 
identify  the  damping  parameter  X4 ,  when  the  initial 
estimate  of  the  parameter  is  in  error  by  100  percent 
are  shown  in  Figure  33  for  two  initial  parameter 
variance  estimates.  When  the  initial  parameter 
variance  is  cased  on  the  initial  estimate  of  X4 ,  the 
near  steady-state  parameter  estimate  is  approximately 
14  percent  high;  for  an  initial  parameter  variance 
that  is  greater  by  approximately  25  percent,  the 
resulting  near  steady-state  value  of  X4  is  nearly 
19  percent  high.  The  corresponding  parameter  variances 
1.5  seconds  after  the  filter  begins  integrating  are 
such  that  the  estimated  values  of  X4  are  within 
approximately  one  standard  deviation  of  the  true  value. 

(4)  Comparison  of  Filter  Performance  to  Chapman-Kirk 
Technique 

The  parameters  of  interest  and  their  uncertainties 
were  extracted  from  the  given  noisy  data  set  using 
the  previously  described  technique  of  Chapman  and 
Kirk  in  order  to  be  compared  to  the  filter  results. 

The  percent  error  in  the  parameters  and  their  normalized 
estimated  standard  deviations  obtained  with  the 
Chapman-Kirk  program  are 


X3:12.5  percent  (1.8  percent) 
X4 :  3.1  percent  (0.3  percent) 


The  above  results  are  similar  to  those  presented  in 
Table  VI  in  that  the  estimated  parameter  uncertain¬ 
ties  are  not  of  sufficient  magnitude  to  reflect  the 
true  error  in  the  extracted  parameters.  Whereas  the 
extended  filter  estimates  the  parameters  and  variances 
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so  that  the  true  value  of  the  parameter  is 
generally  within  one  standard  deviation  of  the 
estimate,  it  is  shown  that  such  is  not  the  case 
when  using  the  Chapman-Kirk  technique.  The  errors 
in  the  above  parameters  are  also  slightly  greater 
than  those  obtained  with  the  extended  filter, 
although  this  is  probably  of  less  consequence  than 
the  differences  in  the  uncertainties  computed  by 
the  two  techniques. 

c.  Nonlinear  System  with  Measurement  Errors  Only 

The  pattern  used  in  the  previous  two  sections  of 
generating  a  data  set,  investigating  the  basic 
response  of  the  filter  in  identifying  the  para¬ 
meters  of  interest,  and  checking  the  sensitivity 
of  the  filter  to  variations  in  the  initial  para¬ 
meter  variances  is  essentially  repeated  here  for 
the  more  complicated  system  model 


a  +  CC4+C6cc2)d  +  (C3+C5C12)  a  =  0  .  (88) 


In  addition,  some  effects  of  errors  in  the  measure¬ 
ment  error  variance  are  also  included  and  discussed. 
Some  comparisons  of  results  achieved  with  the  ex¬ 
tended  filter  and  the  Chapman-Kirk  technique  are 
also  made. 

The  standard  deviation  of  the  measurement  errors  in 
the  generated  data  is  0.00569  radian.  The  mean 
value  of  the  errors  is  zero. 

(1)  Basic  Filter  Performance 

The  time  histories  of  the  four  parameters  of  interest 
as  computed  by  the  extended  filter  while  analyzing 
the  data  just  described  are  shown  in  Figures  34,  35, 
36  and  37.  The  initial  errors  in  the  parameters  are 
-25  percent  for  X3,  X4  ,  and  X6  and  -5  percent  for 
X5.  The  measurement  error  variance  is  that  obtained 
from  the  simulation.  The  initial  values  of  the  para¬ 
meter  variances  are  based  on  the  difference  between 
the  true  values  of  the  parameters  and  the  above- 
mentioned  initial  estimates. 

The  reason  the  initial  estimate  of  X5  is  only  5  per¬ 
cent  low  instead  of  25  percent,  as  is  the  case  with 
the  other  parameters ,  is  due  to  its  magnitude  and  a 
weakness  in  the  numerical  scheme.  The  true  value 
of  X5  is  I960  based  on  the  data  in  Table  III.  A  25 
percent  initial  error  in  this  parameter  results  in 
an  initial  parameter  sample  variance  of 
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p55  (0)  =  (1960-C0. 25)  1960)  2  =  240,000 


which  is  of  sufficient  magnitude  to  cause  the  filter 
to  diverge  for  the  numerical  integration  step  size 
of  0.005  second.  The  divergence  occurs  because  the 
large  initial  value  of  P55  results  in  large  values 
of  P25  and  Pi  5  which,  in  turn,  cause  the  values  of 
P55  that  are  fed  into  the  numerical  integration  sub¬ 
routine  to  be  very  large  [Equations  (81)  and  (82)]. 

The  result  of  this  cascading  effect  is  a  rapid  de¬ 
crease  in  P55  in  large  increments  until  it  becomes 
negative  (a  physical  impossibility  since  this  term 
is  the  variance  of  a  parameter) .  The  structure  of 
the  covariance  equations  is  such  that  a  negative 
variance  causes  numerical  divergence  in  the  solutions 
to  several  of  the  equations . 

To  circumvent  the  problem,  an  initial  error  of 
-10  percent  in  X5,  with  a  corresponding  adjustment 
in  P5 5  (0) ,  was  tried.  The  result  was  again  diver¬ 
gence.  Finally,  an  initial  estimate  of  -5  percent 
was  found  to  be  successful. 

The  results  shown  in  Figures  34,  35,  and  36  indicate 
that  the  filter  does  an  excellent  job  in  correctly 
identifying  the  linear  and  nonlinear  static  restoring 
moment  parameters,  X3  and  X5,  as  well  as  the  linear 
pitch  damping  term,  X4 .  The  results  for  the  non¬ 
linear  damping  parameter,  Xg ,  are  not  as  gratifying 
as  can  be  seen  by  referring  to  Figure  37.  The  error 
in  this  term  is  still  approximately  -10  percent 
after  1.5  seconds.  The  reason  for  this  problem  is 
related  to  the  magnitude  of  this  parameter  as  com¬ 
pared  to  others  in  the  equation  of  motion.  The  true 
value  of  X3  and  X5  are  160  and  I960.,  respectively, 
although  the  magnitude  of  X5  is  effectively  reduced 
by  at  least  an  order  of  magnitude  when  it  is  multiplied 
by  a3.  X4  and  X&,  on  the  other  hand,  have  true  values 
of  1.60  and  4.35,  respectively.  After  X6  is  multiplied 
by  a2d,  it  is  effectively  the  smallest  parameter  in 
the  equation  of  motion  by  an  order  of  magnitude  and 
it  is  naturally  more  difficult  to  identify  since  it 
has  the  least  effect  on  the  trajectory. 

(. )  Effects  of  Variations  in  the  Initial  Parameter 
Variances 

The  effects  of  variations  in  the  initial  parameter 
variances  are  summarized  in  Table  IX  and  in  Figures 
38  and  39.  As  can  be  seen  from  these  results.,  the 
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TABLE  IX.  EFFECTS  OF  VARIATIONS  IN  THE  INITIAL  PARAMETER  VARIANCE 
ON  NEAR  STEADY-STATE  PARAMETERS  AND  THEIR  VARIANCES  (NO 
LINEAR  SYSTEM,  MEASUREMENT  ERRORS  ONLY) 


6] 


relative  to  sample  variance 


near  steady-state  values  of  the  parameters  and 
their  variances  are  relatively  insensitive  over 
the  range  of  initial  parameter  variances  considered. 
Cnly  the  Xg  parameter  estimate  and  its  variance 
show  a  variation  of  more  than  1  percent.  As  is 
obvious  from  the  referenced  table  and  figures ,  no 
data  are  available  for  initial  parameter  variances 
of  50  percent  above  the  sample  variances;  the  reason 
for  this  is,  again,  that  the  magnitude  of  the  P55 
variance  caused  the  filter  to  diverge. 

Attention  is  called  to  the  fact  that,  for  the  range 
of  initial  variances  considered,  the  near  steady- 
state  estimates  of  the  X3,  X4,  and  X5  parameters 
are  all  within  one  standard  deviation  of  their  true 
values  when  the  standard  deviation  is  taken  as  the 
square  root  of  the  near  steady-state  variance  for 
the  parameter  of  interest.  The  Xg  parameter  is 
always  within  two  standard  deviations  of  its  true 
value. 

(3)  Effects  of  Errors  in  the  Estimate  of  the  Measurement 
Error  Variance 

Up  to  this  point,  exact  values  of  the  measurement 
error  variance,  r,  have  been  used  in  all  runs  made 
with  the  extended  filter.  Exact  knowledge  of  this 
quantity  is  available  because  the  pitch-angle  data 
being  used  with  the  extended  filter  are  computer¬ 
generated  with  specified  noise  statistics.  In 
actual  test  situations,  true  values  of  the  noise 
parameters  may  not  be  known  exactly  and,  in  fact, 
methods  for  determining  both  measurement  error 
variances  and  system  noise  variances  from  noisy 
dynamic  data  are  of  current  research  interest. ^4 

Some  effects  of  errors  in  the  estimates  of  the 
measurement  error  variance  on  the  near  steady- 
state  parameter  estimates  and  their  uncertainties 
are  given  in  Table  X  and  in  Figures  40  and  41. 

These  results  show  that  errors  in  the  measurement 
noise  variance  of  ±25  percent  have  essentially  no 
effect  on  estimates  of  the  four  parameters  of  in¬ 
terest.  The  uncertainties  of  the  X3,  X4,  and  X5 
parameters  are  changed  by  less  than  1  percent  over 
the  range  considered.  The  normalized  Xg  uncertainty 
varies  from  approximately  5  percent  to  12  percent. 
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(4)  rilter  Response  to  Larger  Error  in  the  Initial 
Estimate  of  X5 

The  difficulty  associated  with  numerically  large 
initial  parameter  variance  estimates  has  been 
discussed  earlier.  The  sample  variance  of  P55 
which  corresponds  to  a  -5  percent  error  in  the 
initial  X5  parameter  estimate  appears  to  be  an 
approximate  upper  limit  for  the  integration  step 
size  and  physical  constants  used  in  this  analysis. 

A  run  was  made,  however,  with  the  previously 
described  nonlinear  data  set  containing  only 
measurement  errors  where  all  initial  parameter 
estimates  were  in  error  by  -25  percent;  the  P55 
variance,  however,  was  set  equal  to  the  previously 
used  sample  variance  which  corresponds  to  a  -5 
percent  initial  estimate  of  X5 .  The  parameter 
estimate  time  histories  resulting  from  this 
initialization  are  shown  in  Figures  42  through  45. 
These  results  indicate  that  the  near  steady-state 
values  of  the  parameters  are  very  nearly  equal  to 
those  presented  when  the  sample  variances  based  on 
the  initial  parameter  estimates  were  used  initially. 

(5)  Comparison  of  Filter  Performance  to  Chapman- 
Kirk  Technique 

The  linear  and  nonlinear  static  restoring  moment 
and  damping  parameters  were  extracted  from  the 
noisy  data  set  being  considered  in  this  section 
using  the  Chapman -Kirk  program.  The  percent  errors 
in  the  parameters  and  their  normalized  estimated 
uncertainties  are  given  below. 


X3:0.36  percent  (0.35  percent) 
X4t0.69  percent  (0.43  percent) 
Xs:0.41  percent  (0.49  percent) 
Xg:4.60  percent  (11.27  percent) 


Comparing  the  above  results  with  those  previously 
given  in  Table  IX  reveals  that  the  estimates 
obtained  using  the  extended  filter  are  slightly 
more  accurate  for  all  but  the  X4  parameter.  Since 
the  errors  are,  for  the  most  part,  less  than  1  per¬ 
cent,  the  differences  are  essentially  insignificant. 
The  normalized  uncertainties  obtained  with  both 
techniques  are  equally  good. 


d .  Nonlinear  System  with  Measurement  Errors  and 
System  Noise 

In  this  section,  as  before,  the  basic  filter  per¬ 
formance,  sensitivity  of  the  filter  to  initial 
parameter  variance  variations  and  sensitivity  of 
the  filter  to  errors  in  the  estimates  of  the  noise 
variances  are  investigated.  The  data  used  in  the 
analysis  were  generated  by  numerically  integrating 
the  following  nonlinear  equation  with  a  random 
forcing  function, 


a  +  (Cu+C6a2)d  +  (C3+C5ot2)a  =  w(t)  .  (89) 


Approximately  the  same  standard  deviations  used  previously 
are  repeated  once  again;  the  standard  deviation  of  the 
measurement  errors  is  0.00576  radian  and  the  standard 
deviation  of  the  system  noise  is  4.75  rad/sec2. 

(1)  Basic  Filter  Performance 

The  time  histories  of  the  errors  in  the  parameters 
of  interest  obtained  with  the  extended  filter  are 
given  in  Figures  46  through  49.  The  response  time 
required  for  the  parameters  to  reach  their  near 
steady-state  values  is  slightly  greater  than  in  the 
simpler  cases  that  have  been  considered  previously. 
The  accuracy  of  the  parameters  is  similar  to  that 
which  was  achieved  for  the  nonlinear  system  with 
only  measurement  errors  in  the  data  with  the  ex¬ 
ception  of  Xt, ,  which  is  approximately  5  percent 
higher . 

(2)  Effects  of  Variations  in  the  Initial  Parameter 
Variances 

The  effects  of  variations  in  the  initial  parameter 
variances  are  given  in  Table  XI  and  are  shown  in 
Figures  50  and  51.  Only  the  nonlinear  damping  para¬ 
meter  and  its  uncertainty  vary  by  more  than  1  per¬ 
cent  for  initial  variance  variations  of  up  to  50 
percent  relative  to  the  sample  variances . 

(3)  Effects  of  Errors  in  the  Estimates  of  the  Noise 
Variances 

The  effects  of  errors  in  the  estimates  of  r,  the 
measurement  error  variance,  and  q,  the  system  noise 
variance,  are  given  in  Table  XII  and  are  shown  in 
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Figures  52  and  53.  These  figures  show  the 
variation  of  the  near  steady-state  parameter 
errors  and  normalized  parameter  uncertainties 
as  functions  of  the  error  in  the  estimate  oC  the 
measurement  error  variance  for  various  values  of 
error  in  the  system  noise  variance  estimate.  From 
these  results  it  is  obvious  that  only  Xg  is 
affected  by  errors  of  ±25  percent  in  the  estimate 
of  the  system  noise  variance  or  measurement  error 
variance . 

(4-)  Comparison  of  Filter  Performance  to  Chapman-Xirk 
Technique 

The  percent  errors  in  the  parameters  and  their 
normalized  estimated  standard  deviations  that  are 
obtained  using  the  Chapman-Kirk  program  with  the 
data  currently  being  analyzed  are  given  below: 


X3:0.62  percent 
X4:2.63  percent 
X5 :0.41  percent 
Xg : 46 . 3  percent 


(0,40  percent) 
(1.23  percent) 
(0.61  percent) 
(14.8  percent) 


Comparing  these  results  with  the  results  in  Table 
XI  indicates  that  both  techniques  yield  essentially 
the  same  accuracy  in  their  determination  of  X3  and 
X5.  The  filter  does  considerably  better  in  estimating 
Xg  but  is  slightly  worse  in  its  estimate  of  X4 . 

The  uncertainties  computed  using  the  filter  are  of 
sufficient  magnitude  so  that  the  true  value  of  the 
X3 /  X5,  and  Xg  parameters  are  within  one  standard 
deviation  of  the  filter  estimates;  the  true  value 
of  X4  is  within  two  standard  deviations. 

wnen  using  the  Chapman-Kirk  technique,  only  X5  is 
within  one  <=>h-H  mated  standard  deviation  of  its 
true  value;  X3  and  X4  are  within  two  estimated 
standard  deviations  of  their  respective  true  values. 

The  extracted  value  of  Xg  is  not  within  three  estimated 
standard  deviations  of  its  true  value . 


TABLE  XI .  EFFECTS  OF  VARIATIONS  IN  THE  INITIAL  PARAMETER  VARIANCES 
ON  NEAR  STEADY-STATE  PARAMETERS  AND  TKEIR  VARIANCES  'NON 
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.  EFFECTS  OF  LRP.DftS  IN  THE  ESTIMATES  OF  NOISE  VARIANCES 
ON  NEAR  STEADY -STATE  PARAMETERS  AND  THEIR  VARIANCES 
(NONLINEAR  SYSTEM,  MEASDREI !ENT  ERRORS  AND  SYSTEM  NOISE) 


Variation  of  Percent  Error  in  X4  with  Time 
("Linear  System,  Measurement  Errors  Only) 
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Figure  23.  Variation  of  Normalized  Near  Steady-State 

Parameter  Uncertainty  with  Initial  Parameter 
Variance  (Linear  System  Measurement  Errors  Only) 
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Variation  of  X4  Variance  with  Time  for  Three  Initial  Estimates 
(Linear  System,  Measurement  Errors  Only) 


27.  Variation  of  Percent  Error  in  X4  with  Time  (Linear 
System,  Measurement  Errors  and  System  Noise) 
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Figure  28.  Variation  of  Near  Steady-State  Parameter 
Error  with  Initial  Parameter  Variance 
(Linear  System,  Measurement  Errors  and 
System  Noise) 
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Figure  29.  Variation  of  Normalized  Near  Steady-State 

Parameter  Uncertainty  with  Initial  Parameter 
Variance  (Linear  System,  Measurement  Errors 
and  System  Noise) 
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Figure  38.  Variation  of  Near  Steady-State  Parameter  Error 
with  Initial  Parameter  Variance  (Nonlinear 
System,  Measurement  Errors  Only) 
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Figure  39.  Variation  of  Normalized  Near  Steady-State 

Parameter  Uncertainty  with  Initial  Parameter 
Variance  (Nonlinear  System,  Measurement 
Errors  Only) 
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Figure  40. 


Variation  of  Near  Steady-State  Parameter 
Error  with  Error  in  ti e  Estimate  of  the 
Measurement  Error  Variance  (Nonlinear 
System,  Measurement  Errors  Only) 
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!igure  41.  Variation  of  Normalized  Near  Steady-State 
Parameter  Uncertainty  with  Error  in  the 
Estimate  of  Measurement  Error  Variance 
(Nonlinear  System,  Measurement  Errors  Only) 
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Figure  46.  Variation  of  Percent  Error  in  X3  with  Time  (Nonlinear 
Measurement  Errors  and  System  Noise) 
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Figure  50.  Variation  of  Wear  Steady- State  Parameter  Error 
with  Initial  Parameter  Variance  (Nonlinear 
System,  Measurement  Errors  and  System  Noise) 
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Figure  51.  Variation  of  Normalized  Near  Steady-State 

Parameter  Uncertainty  with  Initial  Parameter 
Variance  (Nonlinear  System,  Measurement 
Errors  and  System  Noise) 
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Figure  52.  Variation  of  Near  Steady-State  Parameter 
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Figure  53^  Variation  of  Normalized  Near  Steady-State 
Parameter  Uncertainty  with  Error  in  the 
Estimate  of  Measurement  Error  Variance 
(Nonlinear  System,  Measurement  Errors  and 
System  Noise) 


SECTION  VI 


SUMMARY 


1 .  General  Comments 


An  analysis  of  methods  for  extracting  aerodynamic 
coefficients  from  dynamic  test  data  has  been  conducted. 

The  primary  concern  of  this  analysis  has  been  to  determine 
the  accuracy  with  which  these  coefficients  can  be  determined 
from  dynamic  data  containing  system  noise  and  measurement 
errors.  Both  linear  and  nonlinear  systems  have  been  con¬ 
sidered  with  emphasis  on  the  latter. 

Two  methods  of  extracting  the  coefficients  of  interest 
and  estimating  their  uncertainties  have  been  analyzed.  The 
first  of  these  is  a  deterministic  technique  due  to  Chapman 
and  Kirk. 10  The  second  technique  considered,  based  on  the 
extended  Kalman  filter^ /31  with  augmented  state  vector, 
is  stochastic  in  nature  and  has  been  developed  and  applied 
in  this  report. 

Both  techniques  considered  will  extract  aerodynamic 
coefficients  from  noisy  dynamic  data.  The  degree  to  which 
they  are  successful  has  been  presented  in  detail  in  the 
preceding  chapters. 

2.  The  Chapman-Kirk  Technique 

The  presence  of  noise  in  the  dynamic  data  being  used  by 
the  Chapman-Kirk  technique  affects  the  accuracy  with  which 
this  technique  can  be  used  to  determine  the  coefficients 
of  interest.  When  measurement  errors  only  are  present  in 
the  data,  these  effects  have  been  found  to  be  negligible 
on  linear  and  nonlinear  static  pitching  moment  coefficient 
derivatives;  the  effects  on  the  linear  and  nonlinear  pitch 
damping  coefficients  are  more  pronounced,  especially  for 
the  higher  measurement  errors.  The  error  in  the  linear 
pitch  damping  coefficient  is  as  high  as  25  percent  in  one 
case.  The  greatest  error  in  the  extracted  nonlinear  pitch 
damping  coefficient  is  227  percent.  There  is  an  obvious 
inconsistency  in  the  extracted  values  of  this  nonlinear 
coefficient  since  another  data  set,  which  also  contained 
measurement  errors  with  similar  statistical  properties, 
yielded  a  value  for  this  term  which  was  in  error  by  only 
5  percent. 

The  estimation  of  the  uncertainties  of  the  extracted 
coefficients  when  the  data  contained  only  measurement  errors 
is  quite  adequate.  In  every  case  considered,  the  true  value 
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of  the  coefficient  of  interest  was  found  to  be  within  three 
estimated  standard  deviations  of  the  extracted  value.  The 
estimated  standard  deviations  of  a  given  coefficient  are  also 
consistent  for  different  data  sets  with  the  same  measurement 
error  statistics. 

Use  of  data  containing  measurement  errors  and  system 
noise  has  a  more  adverse  effect  on  the  accuracy  of  coeffici¬ 
ents  determined  with  the  Chapman-Kirk  technique.  As  is  the 
case  when  analyzing  data  containing  only  measurement  errors , 
the  pitch  damping  coefficients  are  the  ones  most  affected. 

The  errors  in  these  coefficients  range  from  1  percent  to 
43  percent  for  the  linear  term  and  from  45  percent  to  600 
percent  for  the  nonlinear  term  over  the  range  of  system  noise 
considered. 

Perhaps  of  even  greater  consequence  than  the  errors  in 
the  damping  parameters  is  the  fact  that  their  estimated 
uncertainties  are  too  small  to  indicate  the  true  errors  in 
the  coefficients  when  the  data  contain  both  system  noise  and 
measurement  errors.  This  is  a  direct  consequence  of  using 
the  root-mean-square-error  between  the  final  solution  to 
the  equation  of  motion  and  the  data  to  calculate  the  parameter 
uncertainties.  This  error  term  is  an  adequate  representation 
of  noise  in  the  data  only  if  all>  errors  are  randomly  super¬ 
imposed  on  the  pitch  angle  output  after  the  equation  of 
motion  has  been  integrated.  The  integrated  effects  of  the 
random  accelerations  representing  system  noise,  however, 
are  such  that  the  final  errors  appearing  in  the  pitch  angle 
due  to  the  system  noise  are  not  necessarily  random  or 
independent  from  one  discrete  time  to  another. 

3 .  Extended  Kalman  Filter 

The  feasibility  of  using  an  extended  Kalman  filter  with 
parameter  augmented  state  vector  for  determining  the  values 
of  aerodynamic  coefficients  and  their  uncertainties  from 
dynamic  test  data  has  been  demonstrated  for  a  one-degree- 
cf -freedom  system. 

The  specific  filter  used  here  generally  reaches  near 
steady-state  estimates  of  the  parameters  in  less  than  one 
second  for  the  system  model  and  error  combinations  considered. 

For  linear  systems  with  measurement  errors  only  in  the  data,  the 
extended  filter  yielded  estimates  of  both  the  static  pitching  moment 
coefficient  derivative  and  the  pitch  damping  coefficient  to  within  1 
percent  of  their  respective  true  values.  Slightly  less  accurate  deter¬ 
minations  of  the  same  two  parameters  are  obtained  using  the  filter  when 
both  measurement  errors  and  system  noise  are  present  in  the  data. 
Variations  in  the  initial  parameter  variances  of  25  percent  seen  to  affect 
essentially  only  the  pitch  damping  term  and  then  only  when  both  measure- 
■f.'ut  errors  and  system  noise  are  present  in  the  data. 
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For  nonlinear  systems  with  measurement  errors  only  in  the  data 
the  estimates  obtained  with  the  extended  filter  for  the  linear  and 
nonlinear  static  pitching  moment  coefficient  derivatives  as  well  as 
the  linear  pitch  damping  coefficient  are  within  approximately  1  per¬ 
cent  or  less  of  their  respective  true  values  for  a  variety  of  initial 
parameter  variances  and  noise  variance  estimates .  The 
error  in  the  estimate  of  the  nonlinear  damping  coefficient 
varies  between  11  percent  and  14  percent.  None  of  the 
coefficients  are  particularly  sensitive  to  variations  in 
the  initial  parameter  variances  or  noise  variances  con¬ 
sidered.  As  is  the  case  for  the  linear  system,  the 
extended  filter  is  less  accurate  in  its  parameter  identi¬ 
fication  when  both  measurement  errors  and  system  noise 
of  the  magnitudes  considered  here  are  present  in  the  data. 
The  linear  damping  coefficient  is  most  affected  by 
having  the  addition  of  system  noise  in  the  data.  The 
error  in  this  coefficient  increases  from  approximately 
1  percent  to  approximately  6  percent  for  the  cases  con¬ 
sidered.  Errors  in  the  noise  variance  estimates  of 
±25  percent  have  little  or  no  appreciable  effect  on  the 
linear  and  nonlinear  static  pitching  moment  coefficient 
derivatives  or  the  linear  pitch  damping  coefficient,  and 
only  a  slight  effect  on  the  nonlinear  pitch  damping 
coefficient. 

For  the  range  of  initial  parameter  errors',  initial 
parameter  variances,  and  noise  variances  considered,  the 
extended  Kalman  filter  produces  an  excellent  estimate  of 
the  parameter  uncertainties  after  integrating  the  Ricatti 
equations  for  approximately  1.5  to  2.0  seconds.  For  every 
case  considered  where  the  initial  parameter  estimates  are 
within  25  percent  of  their  true  values ,  the  near  steady- 
state  parameter  value  is  within  two  standard  deviations 
of  its  true  value,  and  in  the  majority  of  cases  within 
one  standard  deviation,  when  the  standard  deviation  is 
taken  to  be  the  square  root  of  the  near  steady-state 
parameter  variance .  J 


Comparison  of  the  Two  Techniques 


There  is  a  natural  tendency  when  analyzing  two  tech¬ 
niques  to  compare  the  results  achieved  when  solving  similar 
problems.  Both  techniques  investigated  have  strengths  and 
weaknesses . 


From  the  preceding,  discussion  in  this  report,  it  is 
obvious  that  the  Chapman-Kirk  technique  requires  less 
information  initially  than  the  extended  Kalman  filter 
since  it  is  only  necessary  to  make  initial  estimates  of 
the  parameters  in  the  former  whereas  the  parameters ,  their 
variances,  and  the  noise  variances  must  all  be  estimated 
initially  when  using  the  extended  filter. 


There  is  no  provision  to  weigh  the  initial  parameter 
estimates  when  using  the  Chapman-Kirk  technique  whereas 
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this  is  done  through  the  initial  variances  when  using 
the  extended  filter.  This  feature  is  particularly 
important  if  certain  parameters  are  known  with  signifi¬ 
cantly  higher  accuracy  before  the  extraction  process  begins . 
The  Chapman-Kirk  technique  will  sometimes  actually  down¬ 
grade  good  initial  estimates  when  analyzing  noisy  data 
whereas  the  extended  filter  will  tend  to  improve  the 
estimate  or  at  least  not  downgrade  it  for  the  initial 
parameter  variances  considered. 

From  an  economic  point  of  view,  both  techniques  are 
essentially  equal  concerning  computational  costs.  The 
extended  filter  must  numerically  integrate  more  equations 
but  only  has  to  do  this  once;  the  Chapman-Kirk  technique, 
being  an  iterative  scheme,  must  integrate  fewer  equations 
more  times. 

The  programming  requirements  of  the  extended  filter 
are  less  than  those  of  the  Chapman-Kirk  technique  since 
the  former  requires  no  matrix  inversions,  simultaneous 
equation  solutions,  or  logic  for  including  or  excluding 
certain  parameters  from  the  extraction  process . 

5.  Areas  for  Additional  Work 


This  research  is  part  of  a  larger  and  continuing 
program  whose  goal  is  a  critical  study  of  methods  for 
extracting  aerodynamic  coefficients  from  dynamic  data. 

Since  the  research  is  of  a  continuing  nature,  it  seems 
appropriate  to  mention  some  areas  where  additional  work 
should  be  contemplated. 

The  results  presented  here  are  entirely  for  a  one- 
degree-of-freedom  system.  The  effect  of  system  noise  and 
measurement  errors  on  coefficients  extracted  from  coupled, 
multiple-degree-of-freedom  systems  may  prove  to  be 
extremely  interesting  and  is  the  next  logical  area  to  be 
investigated. 

The  response  of  the  extended  filter  to  highly  erroneous 
initial  parameter  estimates  seems  quite  good  for  the  limited 
number  of  cases  investigated.  It  seems  possible  to  take 
advantage  of  this  feature  by  constructing  an  extraction 
algorithm  which  employs  both  an  extended  filter  and  a 
Chapman-Kirk  or  similar  technique  in  tandem  with  the 
extended  filter  providing  the  required  accurate  initial 
estimates  for  the  differential  corrections,  iterative 
technique. 

Additional  work  in  relating  the  effects  of  inadequate 
or  erroneous  modeling  to  the  noise  problem  should  also 
prove  to  be  of  interest. 
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Finally,  when  considering  new  techniques  such  as  the 
use  of  the  extended  Kalman  filter,  it  must  be  concluded 
that  additional  numerical  investigations  will  add  to  the 
confidence  or  more  clearly  define  the  limitations  of  this 
method  of  coefficient  extraction. 
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APPENDIX  I 


CHAPMAN-KIRK  COEFFICIENT  EXTRACTION  NOMENCLATURE 
LIST  AND  PROGRAM  LISTING 


Nomenclature  List 


PROGRAM 

VARIABLE 

H 

ITO 

TMAX 

TZERO 

N 

XZ(1) 


MATH  SYMBOL 
At 


max 


*-0 

N 


a 


o 


DEFINITION 


Numerical  integration 
step  size  (sec) 

Frequency  of  numerical 
integration  output 

Cutoff  time  for  numeri¬ 
cal  integration  (sec) 

Initial  time  for 
numerical  integration 
(sec) 

Number  of  first  order 
differential  equations 
to  be  integrated 

Initial  condition  of  a 
(rad) 


XZ  (2) 


Initial  condition  of  & 
(rad/sec) 


XZ  (3) 


9a 

9Cj 


t=0 


Initial  condition  for 
parametric  differential 
equation 


XZ  (4) 


d_  9g 
dt  9ci 


t=0 


Initial  condition  for 
parametric  differntial 
equation 


XZ  (5) 


9a 

9c2 


t=0 


Initial  condition  for 
parametric  differential 
equation 
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PROGRAM 

VARIABLE 


MATH  SYMBOL 


DEFINITION 


XZ  (6) 


XZ  (7) 


XZ  (8) 


XZ  (9) 


XZ  (10) 


XZ(ll) 


d  3a  Initial  condition  for 

dt  3c2  parametric  differential 

.  n  equation 


Initial  condition  for 
parametric  differential 
equation 


d  3a 
dt  3c3 


Initial  condition  for 
parametric  differential 
equation 


Initial  condition  for 
parametric  differential 
equation 


d  3a 
dt  3c4 


Initial  condition  for 
parametric  differential 
equation 


Initial  condition  for 
parametric  differential 
equation 


XZ (12) 


XZ (13) 


XZ (14) 


XZ (15) 

XZ (16) 


3a 

3c6 

t=0 


d 

dt 


d_  3a 
dt  3C? 

t="0 
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Initial  condition  for 
parametric  differential 
equation 

Initial  condition  for 
parametric  differential 
equation 

Initial  condition  for 
parametric  differential 
equation 

Initial  condition  for 
parametric  differential 
equation 

Initial  condition  for 
parametric  differential 
equation 


I' HOC.  KAN 


VAR]  ABLE 

MATH  SYMBOJ. 

DEFINITION 

XZ (17) 

3  a 

3c8 

t=0 

Initial  condition  for 
parametric  differential 
equation 

XZ (18) 

d  3a 
dt  3ce 

t=0 

Initial  condition  for 
parametric  differential 
equation 

Q 

q 

Dynamic  pressure 
(lb/ft2) 

V 

V 

Freestream  velocity 
(ft/sec) 

AT]  I 

I 

Moment  of  inertia  about 
an  axis  through  the 
vehicle  C.G.  and  normal 
to  the  pitch  plane 
(slug  .  ft2) 

A 

A 

Vehicle  reference  area 
(ft2) 

D 

d 

Vehicle  reference 
length  (ft) 

CMAO 

C 

mao 

Static  pitching  moment 
coefficient  derivative 
(rad-1 ) 

CMA2 

c 

ma2 

Static  pitching  moment 
coefficient  derivative 
(rad-3) 

CMA4 

C 

ma4 

Static  pitching  moment 
coefficient  derivative 
(rad-5) 

CMQO 

Cmqo 

Pitch  damping  coeffi¬ 
cient 

CMQ2 

c 

mq2 

Pitch  damping  coeffi¬ 
cient  (rad-2) 

CMDA 

SnSa 

Static  pitching  moment 
coefficient  at  a  =  0 
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PROGRAM 

VARIABLE 

MATH  SYMBOL 

DEFINITION 

M 

M 

Number  of  experimental 
data  points 

JJJMAX 

Maximum  number  of 
iterations  allowed 
before  program  termi¬ 
nation 

EP 

£ 

Convergence  criteria 
for  change  in  RMS 

ERROR  (rad) 

AS  (I) 

Ci  (t»  ) 

exp 

Experimental  values  of 
a  (rad) 

Cl 

Ci 

a 

o 

C2 

C2 

• 

a 

o 

C3 

c3 

C  q  Ad 

-  mao^1 

I 

C4 

C4 

C  q  Ad2 

mqo 

2VI 

C5 

C5 

C  q  Ad 

-  met? 

I 

C6 

C6 

C  q  Ad2 

"  mq2^ 

2VI 

C7 

c  7 

Cm„  q  Ad 

“  mcti^ 

I 

C8 

c8 

C  .  q  Ad 
moa 

I 


Categories  of  Input  Data 


Numerical  integration  constants,- — These  constants 
include  the  numerical  integration  step  size,  frequency  of 
output,  time  at  which  the  numerical  integration  is  to  stop, 
initial  time,  and  the  number  of  first  order  equations  to 
be  integrated. 

Initial  conditions . — These  data  are  the  first  estimates 
of  the  initial  conditions  for  the  equation  of  motion  as  well 
as  the  initial  conditions  for  the  parametric  differential 
equations . 

Aerodynamic  and  physical  constants. — These  constants 
are  the  dynamic  pressure  and  freestream  velocity  that 
existed  when  the  experimental  pitch  angle  data  were  re¬ 
corded,  the  moment  of  inertia  about  an  axis  normal  to  the 
pitch  plane  and  through  the  center  of  gravity  of  the  vehicle 
of  interest,  a  vehicle  reference  area,  and  a  vehicle  refer¬ 
ence  length. 

Aerodynamic  coefficient  estimates . — These  data  are  the 
first  estimates  of  the  aerodynamic  coefficients  to  be 
extracted . 

Convergency  criteria. — These  constants  are  the  number 
of  experimental  data  points  to  be  read,  the  maximum  number 
of  iterations  to  be  allowed,  and  the  tolerance  against 
which  the  root-mean-square  error  difference  is  tested  for 
convergence . 

Experimental  data. — These  data  are  the  experimental 
pitch  angle  values  from  which  the  aerodynamic  coefficients 
are  to  be  extracted.  The  program  assumes  the  time  between 
each  data  point  is  constant,  and  this  time  must  be  equal 
to  the  product  of  the  numerical  integration  step  size  and 
frequency  of  numerical  integration  output  that  were  given 
on  the  first  data  card. 
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DIMENSION  XZ (18)  ,AE(500)  /X(120)  ,PC (8,500)  ,AC(500)  , 
1AJK (4,4) ,B (8) ,DELC (8) ,L1(4) ,L2(4) , DIF (500) ,STDEV(8) 
2ACDOT (500) 
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WRITE  (6,60)  H . ITO , TilAX , TZERO , N 
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121  WRITE (6, 122) 

122  FORMAT ( 1 0  CONVERGENCE  FAILED  -  MAXIMUM  NUMBER  OF 

1  *  ITERATIONS  EXCEEDED 1 ) 
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1F9.3/5X,  ’CMQO=’  ,F9.3,5X,  ‘CMQ2=’  ,F9.3,5X,  1  CMDA=  1  ,F9.3) 
WRITE  (6,147)  SDCMAO ,  SDCMA2 ,  SDCMA4  ,  SDCMQO  ,  SDCMQ2  ,  SDCMDA 
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SUBROUTINE  XDOT(A,K) 

DIMENSION  A  (18)  ,X(120)  ,PC(8,50C)  ,  AC  (500)  ,ACDOT(500) 
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APPENDIX  II 


EXTENDED  KALMAN  FILTER  NOMENCLATURE  LIST 
AND  PROGRAM  LISTING 


Nomenclature  List 


PROGRAM 

VARIABLE  MATH  SYMBOL  DEFINITION 


H 

TZERO 

TMAX 

ITO 


N 

M 

XZ(I) 

Q 

R 

2(1) 

X(l) 


At 


max 


N 

M 


q 

r 


z 


Xi 


Numerical  integration 
step  size  (sec) 

Initial  time  for  num¬ 
erical  integration  (sec) 

Cutoff  time  for  numerical 
integration  (sec) 

Frequency  of  numerical 
integration  output 

Number  of  first  order 
equations  to  be  integrated 

Number  of  data  points 

Initial  conditions  for 
state  variable  and 
variance  estimates 

System  noise  variance 
(rad2/sec‘1  j 

Measurement  error 
variance  (rad2) 

Experimental  data 
values  of  a 

Estimated  value  of 
a  (rad) 


PROGRAM 

VARIABLE 


MATH  SYMBOL 


DEFINITION 


X(2) 

a 

x2 

Estimated  value  of  d 
(rad/sec) 

X  (3) 

X3 

Estimated  value  of  C3 

X  (4) 

x4 

Estimated  value  of  C4 

X  (5) 

A 

X5 

Estimated  value  of  C5 

X  (6) 

A 

x6 

Estimated  value  of  Cg 

X  (7) 

pll 

A 

Variance  of  Xi 

X  (8) 

p12 

A  A 

Covariance  of  Xi  and  X2 

X  (9) 

P13 

A  A 

Covariance  of  Xj  and  X3 

X(10) 

p14 

A  a 

Covariance  of  Xi  and  X4 

x(ll) 

p15 

A  A 

Covariance  of  Xj  and  X5 

X  (12) 

P1  6 

A  A 

Covariance  of  Xj  and  X6 

X  (13) 

P22 

A 

Variance  of  X2 

X  (14) 

P23 

A  A 

Covariance  of  Xj  and  X3 

X  (15) 

p  2  4 

A  A 

Covariance  of  Xi  and  X4 

X  (16) 

p  2  5 

A  A 

Covariance  of  X2  and  X5 

X  (17) 

P2  6 

•  A  A 

Covariance  of  X2  and  Xg 

X  '  1 8 ) 

P  3  3 

A 

Variance  of  X3 

X  (19) 

P  3  4 

A  A 

Covariance  of  X3  and  X4 

X(20) 

P  3  5 

A  A 

Covariance  of  X3  and  X5 
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PROGRAM 

VARIABLE 

MATH  SYMBOL 

DEFINITION 

X  (21) 

p36 

■a  a 

Covariance  of  X3  and  X6 

X  (22) 

P44 

A 

Variance  of  X4 

X  (23) 

p45 

A  A 

Covariance  of  X4  and  X5 

X  (24) 

p46 

A  A 

Covariance  of  X4  and  Xg 

X  (25) 

p55 

A 

Variance  of  Xs 

X  (26) 

P56 

A  A 

Covariance  of  X5  and  X6 

X  (27) 

P6  6 

Variance  of  X6 

A  (I) 

• 

X 

Derivative  of  X  with 
respect  to  time 
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Categories  of  Input  Data 


Numerical  integration  constants  ,—^These  constants  are 
the  numerical  integration  step  size,  frequency  of  output, 
time  at  which  the  integration  is  to  stop,  initial  time,  the 
number  of  first  order  equations  being  integrated,  and  the 
number  of  data  points  to  be  read. 

Initial  conditions. — These  data  are  the  initial  estimates 
of  the  state  variable  components ,  state  variable  variances ,  and 
covariances . 

Noise  variance  estimates. — These  two  entries  are 
the  estimates  of  the  measurement  error  variance  and  system 
noise  variance. 

Pitch  angle  data. — These  data  are  the  experimental 
values  of  the  pitch  angle  from  which  the  filter  estimates 
the  states  and  parameters.  The  time  between  data  points 
is  not  input  but  is  assumed  to  be  the  same  as  the  numerical 
integration  step  size  input  on  the  first  data  card. 
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